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ABSTRACT 


The  experimental  study  of  the  radiative  capture  of  protons  by 
tritium  was  performed  at  proton  energies  of  227  MeV,  300  MeV  and 
375  MeV.  A  polarized  beam  was  used  and  measurements  were  done  at  a 
number  of  angles.  The  cross  section  and  the  analyzing  power  of  the 
reaction  were  obtained. 

The  study  of  the  reaction  beyond  meson  production  threshold 
enabled  the  Meson  Exchange  effects  to  be  seen  explicitly.  These  ME 
effects  are  implicitly  included  in  the  low  energy  region  through 
electric  dipole  operators  by  means  of  the  Siegert ’s  Theorem.  First 
order  correction  through  Gauge  invariance  gives  rise  to  the  seagull 
term  which  shows  up  as  a  NNtty  four-point  vertex  term.  The  other 
factors  contributing  to  the  cross-section  are  the  A(1232)  and  short- 
range  correlations  generating  from  the  high  momentum  transfer 
characteristic  of  the  (p , y)  process. 

Phenomenological  models  like  the  Distorted  Wave  Impulse 
Approximation  (DWIA)  are  quite  successful  in  describing  the  (p ,  y) 
process.  DWIA  is  an  extension  of  the  Quasideuteron  model  introduced 
by  Levinger.  The  interaction  vertex  is  replaced  by  the 
experimentally  known  cross  section  of  p+n  ->  dy.  The  additional 
factors  coming  into  the  expression  for  the  cross  section  are  a  form 
factor  which  is  the  fourier  transform  of  the  overlap  of  the 
wavef unctions  with  respect  to  the  momentum  transfer,  a  spin  factor 
and  a  kinematical  factor. 
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Another  equally  important  motivation  in  studying  the  reaction 
was  to  check  the  validity  of  the  Time  Reversal  Invariance  (TRI)  in 
the  electromagnetic  interaction  of  the  hadrons.  The  principle  of 
detailed  balance  relates  the  cross  section  for  the  forward  reaction 
to  that  of  the  inverse  reaction  on  the  assumption  of  TRI. 

Comparison  of  the  cross  section  values  obtained  in  the  present 
experiment  with  the  inverse  reaction  data  provided  a  check  for  the 
validity  of  TRI. 

Comparison  of  the  results  with  the  theoretical  predictions  by 
Gari  and  Heback  revealed  the  significance  of  the  A(1232) 
contribution.  The  DWIA  calculation  agrees  with  the  results  very 
well  except  at  high  energies  and  backward  angles.  This  points 
towards  the  multi-step  processes  which  become  important  at  high 
energies  and  backward  angles.  Comparison  of  the  data  with  the 
inverse  reaction  data  reveals  very  good  agreement  between  them,  both 
in  the  magnitude  of  the  cross  sections  and  the  shape  of  the  angular 
distributions,  confirming  the  validity  of  TRI  in  the  electromagnetic 
interaction  of  the  hadrons. 
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CHAPTER  1 


INTRODUCTION 

This  thesis  reports  an  experimental  study  of  the  reaction 
P  +  T  -►  4He  +  y  at  proton  energies  of  227  MeV,  300  MeV  and  375  MeV. 
The  experiment  was  done  at  a  number  of  scattering  angles  to  study 
the  angular  distribution  of  the  cross-section.  Polarized  proton 
beam  was  used,  thus  providing  us  with  the  analyzing  power  of  the 
reaction  as  well. 

The  study  forms  part  of  a  larger  program  of  using  the 
electromagnetic  probe  to  look  at  light  nuclei  such  as  %  and 

4He.  This  has  been  motivated  by  the  desire  to  gain  information  on 
the  nuclear  structure  as  well  as  the  reaction  mechanism. 

Furthermore,  comparison  of  the  reaction  with  the  existing  inverse 
reaction  data  (Pic-70,  Kie-73,  Arg-75,  Are-70,  Sch-83)  allows  a  test 
of  one  of  the  fundamental  symmetries  of  nature,  Time  Reversal 
Invariance  (TRI),  in  electromagnetic  interaction  of  hadrons. 

1.  1  THE  SIGNIFICANCE  OF  THE  REACTION 

The  electromagnetic  nature  of  the  interaction  provided  an 
excellent  tool  to  probe  the  nucleus.  As  it  is  a  weak  interaction 
compared  to  the  nuclear  strong  interaction,  it  acts  as  a 
perturbation  to  the  nuclear  system.  Moreover,  the  electromagnetic 
interaction,  unlike  the  strong  interaction,  is  well  understood 
through  quantum  electrodynamics.  The  (p,y)  reaction  is 
characterized  by  high  momentum  transfer  to  the  residual  nucleus. 

This  means  high  momentum  components  of  the  nuclear  wavefunction  are 
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sampled.  We  get  an  opportunity  to  look  at  the  exotic  degrees  of 
freedom  of  the  nuclear  wavef unction  like  short-range  correlations. 

As  we  go  up  in  the  energy  scale,  various  reaction  mechanisms  make 
themselves  visible  in  varying  orders  of  significance.  The  Direct 
knockout  term  also  known  as  the  Born  terra  or  Direct  reaction  term, 
dominant  at  lower  energies,  becomes  less  important  at  high  energies. 
As  the  energy  is  increased,  the  Two  Nucleon  Model  (TNM)  term,  also 
known  as  the  two  step  mechanism,  phenomenologically  described  under 
Quasideuteron  model,  becomes  important.  Even  though  at  lower 
energies  this  effect  is  buried  in  the  electric  dipole  transition 
through  Siegert’s  theorem,  as  the  energy  increases  it  becomes 
necessary  to  include  it  explicitly  through  meson  exchange  currents. 
Mesonic  and  isobaric  degrees  of  freedom  become  apparent.  As 
described  in  Gar-81,  the  dominant  mechanisms  are  described  by  the 
Shell-Model  terra,  the  Meson  Exchange  Current  (MEC)  term,  the 
Correlation  terra,  and  the  A  isobar  production  term.  In  addition  to 
looking  at  various  reaction  mechanisms,  we  also  get  an  opportunity 
to  look  at  the  nuclear  structure  aspects  of  the  nuclei  %  and  *fae. 

Since  the  reaction  is  related  to  y  +  ^e  P  +  T  through  the 
principle  of  detailed  balance,  we  need  to  discuss  only  one  of  the 
reactions.  In  the  following  discussions  we  will  be  concentrating  on 
the  reaction  y  +  ^e  -►  P  +  T. 

1.2  VARIATION  OF  PHOTODISINTEGRATION  CROSS  SECTION  AS  A  FUNCTION  OF 

ENERGY 

The  total  cross  section  for  the  photodisintegration  is 
characterized  by  a  peak  at  the  incident  y  energy  between  the 
threshold  and  40  MeV.  This  peak  is  due  to  giant  resonance.  The 
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integrated  absorption  cross  section  in  this  peak  is  a  large  fraction 
of  the  value  expected  from  the  total  nuclear  dipole  oscillator 
strength.  Beyond  the  giant  resonance  region,  the  cross  section 
falls  off  monotonically .  Around  300  MeV,  the  cross  section 
increases  due  to  the  A  isobar.  This  is  a  nucleon  excitation  with  a 
free  width  of  about  115  MeV.  Meson  exchange  currents  play  an 
important  role  in  this  energy  region. 

The  low  energy  cross  section  may  be  expressed  in  terms  of 
multipole  contributions,  without  being  concerned  with  the  meson 
exchange  current  or  isobar  production.  The  total  Hamiltonian  is 
expressed  as  the  sum  of  contributions  from  the  nucleus,  from  the 
radiation  field  and  from  the  interaction  between  them.  Using  the 
time-dependent  perturbation  theory,  the  probability  of  transition 
from  the  eigenstates  of  the  uncoupled  system  to  the  eigenstates  of 
the  coupled  system  is  found.  For  a  quantum  system  the  interaction 
of  the  radiation  is  given  as  the  scalar  product  of  two  operators, 
one  refering  to  the  radiation  field  and  the  other  to  the  nucleus. 

The  Hamiltonian  of  the  radiation  field  is  given  by: 


H 

Y 


Z 

kX 


\  alx\\ 


(1.1) 


where 


akXand  akX 


are  the  creation  and  annihilation  operators  for 


photons  with  wavevector  k  and  circular  polarization  X  =  (±1), 
respectively,  a^ ^  a^^  gives  the  number  of  photons  of  type  (k,  X) 


The  radiation  field  is  characterized  by  a  vector  potential  A  given 


by: 
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*  „  r  2-nhc  f  +  ilt-r  ,  ,  .  /io\ 

„ ST"  (akAekxe  +  h'c*)  (u2) 

k,  A  ^ 

where  V  is  the  normalization  volume  and  e^^  the  polarization 
vector. 

The  interaction  Hamiltonian  is  given  in  terms  of  the  vector 
potential  A,  nuclear  current  density  j”  and  nuclear  magnetization 
in  the  following  fashion. 


Hi=-/(^j>.A+J*^xA)d3r  (1.3) 

The  total  Hamiltonian  for  the  electromagnetic  interaction  is 
represented  by  H: 


H  =  H  +  H  +  H.  (1.4) 

n  r  l 

where  H^  represents  the  Hamiltonian  of  the  nucleus.  The  unperturbed 
system  Hamiltonian  is  represented  by  H^  which  equals  H^  +  H^.  Let 
the  uncoupled  system  eigenstates  be  represented  by  |N>.  The 
interaction  H^  is  turned  on  at  t  =  to  and  the  eigenstates  | N> 
develop  into  eigenstates  |  4^(t)>  of  the  total  system.  This  is 
expressed  in  the  following  expansion  of  4^  along  the  basis  |N’>  with 
time  dependent  expansion  coefficients  aN,N(t): 

V 

IVt)>=  1  lN’>aN'N(t)  exp 
N 1 


(1.5) 


5 


The  probability  per  unit  time,  WN»N  for  the  transition  from  N 
to  N’  is  expressed  in  terms  of  a^f^  and  the  density  of  final  states 
P(EN,): 

Vn  -5F  /  laN'N(t)l2  p(EN->  dV  0.6) 


Using  eqn.  (1.5)  in  the  time-dependent  Schro'dinger  equation, 
we  solve  for  a^,^.  This  is  substituted  in  eqn.  (1.6)  to  yield  the 
transition  probability: 


NfN 


2  7T 
■ft 


-  <N ' |H± |N>  + 


fi 


s  ir~  <N’  |H  |N"> 
N"  N"N 


<N” |H± | N> 


(En.) 


(1.7) 


The  interaction  Hamiltonian  may  be  expressed  in  terms  of  the 
angular  momentum  transferred  to  the  nucleus.  The  nucleus  may  be 
described  in  terms  of  the  projection  of  nuclear  spin  on  the 
quantization  axis  (quantization  axis  standing  for  the  propagation 
vector  of  the  incident  photon  beam). 


H£  =  e  /  ¥*  zl  »  (1.8) 

e  i  e 

H£  =  e  /  y*  on  V  d  ^r  (1.9) 

m  J  i  £  e 


where 

e 


is  the  matrix  element  for  the  transition  between  the 
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initial  state,  Y  ,  and  the  final  state,  4^.  The  suffix  l  stands  for 
the  electric  multipole  transition  i.  The  matrix  element  for  the 
magnetic  multiple  transition  is  H^. 

The  first  extensive  calculation  for  the  transition  probability 
in  terms  of  various  multipoles  was  done  by  Partovi  (Par-64)  in  the 
context  of  photodisintegration  of  deuterium,  Partovi  applied 
nonrelativistic  phenomenological  Schro'dinger  theory  to  calculate 
the  transition  probability.  Approximations  made  are  the  neglect  of 
nucleon  structure,  pion  exchange  currents  and  multipoles  higher  than 
the  octupole.  The  nucleons  were  assumed  to  interact  through  a 
semiphenomenological  potential  developed  by  Hamada  and  Johnston 
(Ham-62)  given  in  the  following  form: 

V  -  Vc  +  VXS12  +  VLS  (L+  •  S+)  +  VLL  L12  t1*10) 

where  V  is  the  central  potential,  V  the  tensor  potential,  V  the 

C  1  tio 

spin-orbit  potential,  and  V  the  quadratic  spin  orbit  potential. 

The  vector  potential  of  the  radiation  field  was  expanded  in  plane 
waves  of  circular  polarization. 

A(x)  =  Z  Z 

oj  y=±l 

where  a  is  the  annihilation  operator  for  a  photon  of  momentum  a) 

% 

and  polarization  y.  The  interaction  Hamiltonian  is  H?. 

H’  =  -  /  J  •  A  d3r 


(2  tt  i  ^  22 
'  0)  ' 


+  iw*x  T  •> 
a  e  e  -  a  e  e 
w  y  03  -y 

y  y 


->•  -► 
■i  w*x 


(l.n) 


(1.12) 
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where  J  is  the  current  density  operator  of  the  n-p  system.  The 

g 

final  state  is  characterized  by  total  spin  S  and  projection  m  of 
the  n-p  system  and  the  initial  state  by  the  photon  polarization  y 
and  the  deuteron  projection  quantum  number  ra^.  The  transition 
probability  between  these  states  was  found  using  the  interaction 
Hamiltonian  given  by  eqn.  (1.12). 

Partovi's  calculation  covered  energy  regions  up  to  140  MeV.  In 
this  region  the  photodisintegration  cross  section  may  be  fitted  with 
a  polynomial  of  the  form  given  below: 

a  sin20  (1  +  b  cos  0  +  c  cos  20  +  d  cos  30  +  e  cos  40)  (1.13) 

d  u 

where  0  is  the  proton  cm  angle.  Eqn.  (1.13)  is  an  electric 
multipole  expansion  including  El,  E2  and  E3  transitions.  The 
electric  dipole  term  is  dominated  by  the  transition  %  -  3P. 

1.3  MESON  EXCHANGE  THEORIES 

The  experimental  results  of  higher  energies  could  not  be 
explained  if  multipole  transitions  alone  were  taken  into  account. 
Explicit  inclusion  of  meson  exchange  current  along  with  the  one-body 
current  (or  convection  current)  is  necessary  to  take  care  of  the 
total  nuclear  current.  The  total  current  thus  obtained  modified  the 
interaction  Hamiltonian  which  yields  additional  terms  in  the 
transition  matrix.  Furthermore,  imposition  of  the  gauge  invariance 
condition  yields  additional  terms  in  the  transition  matrix.  In  the 
low  energy  situation,  meson  exchange  effects  are  taken  into  account 


to  first  order  implicitly  through  an  electric  dipole  term  applying 
current  conservation.  This  is  done  by  Siegert *s  theorem.  We  first 
take  a  look  at  Sieger t’s  theorem. 
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The  total  nucleon  current  ^  is  given  in  terms  of  the  convection 
current  and  the  meson  current  J^. 


•f- 


(1.13a) 


This  is  substituted  in  eqn.  (1.12)  to  get  the  interaction 

■> 

Hamiltonian.  The  vector  potential  A  is  substituted  in  terras  of  the 
polarization  vector  e  and  the  wave  vector  k.  We  get: 


-►  -»• 

tt  r  t  *  ik  *r  ,  -► 

H.  =  J  •  £  e  dr 
mt  J 


(1.14) 


•>  ik  *r  rj  /  +  +  ik*rv  f>  ,  +  ->v  ik  «r 
ee  =  V  (  e*r  e  )  -  ik  ( e*r)  e 


(1.15) 


In  the  long  wavelength  limit,  the  second  term  in  eqn.  (1.15) 
vanishes  and  we  get: 


H.  ~  /  J  •  V  (e*r  eik  r)  dr^ 
int  J 


(1.16) 


Performing  the  integration  by  parts  and  assuming  *f(r)  goes  to 
zero  for  large  r,  we  get: 


Hint~-/(V.J+)(^)eikVdr 


(1.17) 
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Now  we  invoke  the  equation  of  current  conservation: 

-  0  0.18) 

where  p  is  the  charge  density  contributed  both  by  the  nucleon  field 
and  the  meson  field.  Rewriting  the  above  equation,  we  get: 


H. 

int 


~  / 


3p 


(  e-r) 


ik 


dr 


(1.19) 


Now  making  use  of  Heisenberg's  equation  of  motion: 


3p 

3t 


=  £  [H,p] 


(1.20) 


Using  the  above  two  equations  and  the  fact  that  the  initial 
state  | i>  and  the  final  state  | f >  are  eigenstates  of  H  with 
eigenvalues  and  E^,  we  get: 


<flHint|i>  (W  <f|p0^r+  e  1°  (1-21) 

where  pQ  is  the  expectation  value  of  p.  Comparing  eqns.  (1.14)  and 
(1.21),  we  see  that: 


•f  •  e  =  ^  (Ef-E1)po  t  •  r  (1.22) 

This  means  that  in  the  long  wavelength  limit,  the  matrix 
element  of  the  interaction  Hamiltonian  is  proportional  to  the  matrix 
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element  of  the  electric  dipole  moment  p^r.  This  is  known  as 
Siegert  's  theorem.  Note  that  this  applies  in  the  long  wavelength 
limit  or  low  energy  region. 

Because  Siegert’s  theorem  is  valid  to  first  order  only,  it 
cannot  take  care  of  the  reaction  mechanism  at  higher  energies. 
Explicit  inclusion  of  the  resonant  and  non-resonant  meson  exchange 
current  terms  is  necessitated  at  higher  energies.  Calculations  in 
this  regard  have  been  done  (Lon-76,  Fin-76,  Lon-79,  Gar-81, 

Aren-82).  Described  below  are  the  salient  features  of  the  MEC 
theory. 

We  start  with  the  time-dependent  Schro’dinger  equation  and  see 
how  the  gauge  transformation  for  the  vector  potential,  scalar 
potential  and  the  wavefunction  yields  the  perturbation  expansion  for 
the  interaction  Hamiltonian. 

-  ["IE  ($-^it)2  +  e*  +  v]  f  (1.23) 

the  above  is  the  time-dependent  Schro*dinger  equation  in  the 

presence  of  the  electromagnetic  field  characterized  by  the  vector 
-► 

potential  A  and  scalar  potential  <j>.  Now  we  apply  gauge 
transformation  as  follows: 

A  +  k  +  $X  (1-24) 


X 


(1.25) 
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(1.26) 


Substitution  of  eqns.  (1.24),  (1.25)  and  (1.26)  into  (1.23) 
should  retain  the  form  of  the  Schrodinger  equation  intact  in  order 
for  gauge  invariance  to  hold  good.  This  yields  the  following 
equation  for  the  Hamiltonian: 


H(  C.x) 


^  x 
e 


H(0 


le 


(1.27) 


where  H  is  the  Hamiltonain  and  £  are  internal  operators.  We  see  the 
RHS  of  eqn.  (1.27)  may  be  expanded  in  the  following  fashion: 

H(C,x)  =  H(0  +  H'(c,x)  +  ...  (1.28) 

where  H  is  the  unperturbed  Hamiltonian  and  H'  etc.  are  corrections 
imposed  by  the  gauge  invariance.  We  may  also  observe  the  following 
equation  from  eqn.  (1.27): 

H’(C,X)  =  [  X>H]  (1.29) 

The  total  nuclear  current  is  given  as  a  sum  of  one-body  current 
and  exchange  current: 


J  +  J, 


(1.30) 
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where : 


J  =  E 
o 


[ 


] 


(1.31) 


which  is  the  one-body  current.  The  first  term  refers  to  the 
convection  current  and  the  second  terra  to  the  magnetization  or  spin 
current.  The  meson  exchange  current  J„  stands  for  two,  three,  etc. 

£!i 

body  currents.  In  the  approach  by  Gari  and  Heback  (Gar-81)  the 
nuclear  Hamiltonian  H  is  split  into  a  one-body  operator  T  and  two- 
body  operator  V. 


(1.32) 


H  =  T  +  V 


The  requirement  of  current  conservation  for  the  total  current  JT 


yields  two  relationships,  one  for  the  one-body  current  and  the 
other  for  the  exchange  current  J„. 

Ci 


*  •  J"o  -  'S  [T>p] 


(1.33) 


*  •  =  “S'  Iy*P] 


(1.34) 


where  p  is  the  charge  density.  The  form  of  H  in  eqn.  (1.32)  is 
replaced  by  another  form  as  shown  below. 


H  =  (T  +  U)  +  (V  -  U)  =  H' 


.SM 


+  R 


(1.35) 


SM 

where  U  is  the  single  particle  potential,  H  the  shell-model 
Hamiltonian  and  R  the  residual  interaction.  U  is  taken  to  be  of  the 
Wood-Saxon  type: 


U 


-U 

_ o _ 

1  +  exp  ((r-rQ)/a  J 


(1.36) 


where  U  =58.5  MeV,  r  =  2.77  fm  and  a  =  0.5  fm  for  ^e. 
o  o 

The  nuclear  states  |  y>  are  expanded  in  terms  of  shell-^nodel 
eigenstates  |  <J>>. 


+ 

X"> 


f 


R 

H  +  ie 


(1.37) 


*i>  = 


>  + 


PR 


E. 

l 


rSM 


k> 


(1.38) 


where: 


P  =  1  -  |  ^><*.1  (1.39) 

Now  it  is  possible  to  express  the  transition  matrix  between  the 
nuclear  states  in  terms  of  shell-model  states.  The  transition 
matrix  is  represented  by  M. 


Mfi  =  <Xf  I  “  c  /  ^  (*)  *  •tp  (*)  d3x|  xi> 


(1.40) 


Substituting  eqns.  (1.37),  (1.38)  and  (1.39)  in  the  above 


equation,  Gari-Heback  find  that: 


fi 


=  < 


<J>f  I  ^1  $± 


>  +  < 


Ef  -H 


+  i  e 


a|  <f>  >  +  <  <}>  |  8 


PR 


Ei  -H 


SmI  Xi> 


+  <<t>f 


R 

Ef  -  H  +  ie 


a 


E. 

i 


PR _ 

-hsm 


♦<> 


(1.41) 


where: 


a  =  -  ~  /  A  (x)  •  J_  (x)  d  ^ 


a  =  a  +  a 

o  exc 


(1.42) 


where  a^  is  the  one-body  contribution  and  ^exc  is  the  exchange 
contribution. 

aQ  is  contributed  by  the  one-body  convection  current  and  one- 
body  spin  current  as  shown  below: 


a  = 

o 


2-nfric  ^ 
k 

Y  OF 


>  a  it  >  ,  -►  >  x 

e  e  •  — - -^-r  p  a  •  (  e.x  k  ) 

L  a  A  M  2M  a  a  A  y  J 


ik  «r 
Y  a 


(1.43) 


Let  us  take  a  look  at  each  terra  in  (1.41).  The  first  terra  is 
split  into  two  terras  through  eqn.  (1.42),  one  of  them  standing  for 
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the  shell-model  term  and  the  other  one  standing  for  the  exchange 
current  term.  Note  that  in  the  long  wavelength  limit  the  exchange 
current  term  becomes  the  Siegert  term  we  saw  earlier. 

Gexc  =  [V,Q]  (1.44) 

where  Q  is  the  sum  of  all  the  electric  multipole  operators  Q^. 


Q  =  l 

L=1 


“  2  Tike  f2L+l  >!1/2 

L=1  /IT 
Y 


+  r 


i.  > 

dr  ' 


(>Y> 


ylx  (r) 


] 


(1.45) 


The  two-body  potential  V  related  to  the  exchange  current 
through  current  conservation  is  given  by  the  following  expression: 


“V 

v  =  -  Vo  ^ —  [a  0  +  ajOj.02  +  aTT1«T2  +  a  aT(  ’<J2)  (  T1  ’xi)  1  d*46) 


where  V  =95  MeV,  an  =  a  =  -0.0025,  a  =  -0.1025,  and  a 
-0.2325,  for  4He. 

The  second  and  third  terms  in  eqn.  (1.41)  are  known  as 
correlation  terms,  arising  from  the  residual  interaction.  These  are 
described  by  giant  resonances  of  L=1  and  L=2  for  **He.  The  last  term 
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does  not  contribute  in  the  first  order. 

The  A(1232)  resonance  was  introduced  as  part  of  two-body 
current : 


Res  -y  r 

A  ^ 

A _ \  e 

pRes  J  u 

y"ea 

The  terms  mentioned  so  far  are  pictorially  represented  as  in 
fig.  1.1a  and  fig.  1.1b.  Fig.  1.1a  gives  the  one-body  interaction 
term  and  fig.  1.1b  gives  the  exchange  current  terms.  The  first  term 
in  the  exchange  current  contribution  is  the  meson  exchange  current 
(MEC)  term,  the  second  the  seagull  terra,  and  the  third  the  A  isobar 
term. 

Gari-Heback  found  that  with  ^He  there  is  a  considerable 
contribution  from  exchange  currents  in  the  energy  region  40-400  MeV. 
The  correlation  contribution  was  found  to  decrease  as  energy 
increased.  The  shell  model  term  showed  no  structure  at  30°  whereas 
it  has  a  minimum  at  120  MeV  for  90°. 

The  A  isobar  contribution  was  calculated  in  detail  by  Londergan 
et  al.  (Lon-76)  in  the  general  case  of  a  nucleus  with  A  nucleons. 
They  considered  the  diagram  in  fig.  1.2  for  A  formation. 

The  S  matrix  for  the  reaction  was  given  by: 

S  =  (2*)**  6-  (kA  +  ky  -  kN  -  kA_p  M 


(■^xip 3  (k  x  (apcSp) 


(1.47) 


0.48) 
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FIGURE  1.1a 

ONE-BODY  INTERACTION  DIAGRAMS 


FIGURE  1.1b 


MESON  EXCHANGE  DIAGRAMS 
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A  ISOBAR  DIAGRAM 


where  k's  are  various  momenta  and  M  the  amplitude  for  transition,  M 


was  calculated  by  taking  into  account  the  contributions  at  the 
vertices  ttNA,  ttNN  and  AyN.  The  following  expression  was  obtained 
for  the  transition  amplitude. 

<A-1|h.,m  I  (A-l)*irXP tt| H  I  A> 
_ 1  NNtt1  _ 1  NAtt1 _ 

[%  A“p y)~  wA~e  Y-i  ^  [Ea^A^^A”1^ 

states 

<AIH-YNaIN’Y><N(A_1)*IHnIA> 

[E(A-1)*  (ky-PA>  +  “A  -  E  -  -  “A  -  EN  (V*  ' 

(1.49) 

where  is  the  y  energy,  and  E(P)  is  the  total  energy  corresponding 

* 

to  the  momentum  P.  The  term  (A-l)  refers  to  the  intermediate 
state.  The  nuclear  wavefunctions  were  taken  as  the  harmonic 
oscillator  ones.  Photon-nucleon  vertex  was  calculated  using  both 
convection  and  magnetic  moment  contributions.  The  NAy  vertex  was 
computed  using  the  quark  model  of  N  and  A  as  well  as  QED.  The  it 
vertex  was  described  by  the  spin-3/2  Rarita-Schwinger  field,  Hf 
the  nucleon  field,  and  <f>  the  pion  field. 

* 

Haxt  =  -i  (— )  fT  f  3  (b  (1.30) 

ANtt  m  \i  y 

TT 

•k 

where  T  is  the  baryon  isospin  operator  and  f  =  2.097  determined 


M  =  I 

intermediate 


dP  Ad  0). 
A  A 

i(2 tt)  4 


from  A  decay  width.  Pion  absorption  was  taken  into  account  by 
summing  NNtt  vertices  over  the  nucleons  in  (A-l)  and  A: 
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(1.51) 


71 


where  f  =  1.009. 

The  form  factors  were  included  at  the  ttNN  and  ttN A  vertices  as 
shown  below: 


(1.52) 


where  X  =6. 8m,  A...  =  6. 9m. 

mT  TT  AN  TT  77 

NNtt 

More  recently  Arenhovel  (Are-82)  calculated  the  contribution  of 

meson  currents  and  isobar  propagation.  Exchange  current 

contributions  to  various  electromagnetic  processes  in  nuclei,  mostly 

to  magnetic  transition  Ml,  were  studied. 

The  inclusion  of  N'(1470)  in  the  ground  state  of  **He  was  found 

to  be  useful  in  getting  the  photodisintegration  cross  section  at 

<  400  MeV  (Lag-75).  This  study  was  motivated  by  the  fact  that 

even  though  A(1232)  photoexcitation  is  the  dominant  process  in  the 

reactions  d(y,p)n  and  4He(y,p)3H,  the  cross  section  for  the  former 

is  about  40  times  that  for  the  latter  for  E  =  280  MeV  at  0  =  90°. 

Y  P 

The  cross  section  was  expressed  in  the  following  fashion: 
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d  a 
3TI  a 


/l-  3 2  (A+B)  +  3C 


2 


where : 


%e  =  /1_62  ^  +  0 

where  $2  is  the  probability  of  finding  a  Nf(1470)  in  the  ‘He 
nucleus*  Terms  A,  B  and  C  are  represented  by  the  diagram  in 
fig.  1.3  Laget  found  that  a  3%  admixture  of  N'(1470)  in  the  ground 
state  of  4He  reproduced  the  cross  section  at  medium  angles. 

Arends  (Are-78)  notes  that  contrary  to  photodisintegration  of 
the  deuteron,  but  similar  to  the  reaction  ^eCy,?)^!,  there  is  no 
obvious  evidence  for  a  A(1236)  resonance  contribution  to  the 
reaction  t+He(y,p)%.  Theoretical  investigation  of  the  resonant  and 
non-resonant  pion  exchange  effects  was  done  by  Finjord  (Fin-76). 
Calculations  were  done  for  the  energy  region  165  MeV  to  330  MeV. 

Only  the  s-parts  of  the  wavef unctions  were  considered  and  the  final- 
state  interactions  between  the  proton  and  nucleus  were  neglected. 
Finjord  found  that  non-resonant  pion  exchange  contributions  are 
dominant.  However  his  results  are  strongly  dependent  upon  the 
asymptotic  form  of  the  momentum  space  wavef unction. 

1.4  QUASIDEUTERON  AND  OTHER  PHENOMENOLGICAL  MODELS 


Levinger  (Lev-51)  proposed  the  quasideuteron  model  to  explain 


> 


FIGURE  1.3 


PROTON  POLE,  TRITON  POLE,  Nf(1470)  DIAGRAMS 


TERM  A  :  PROTON  POLE 


TERM  B  :  TRITON  POLE 


TERM  C  :  N’ (1470) 


the  cross  section  of  nuclei  with  A  >  3.  He  expressed  the  nuclear 
wavefunction  in  terms  of  the  deuteron  wavefunction  on  the  assumption 
that  two  very  close  nucleons  inside  the  nucleus  constitute  a 
quasideuteron.  The  wavefunction  of  the  nucleus  was  expressed  as  a 
product  of  the  deuteron  wavefunction  and  the  wavefunction  for  the 
remaining  nucleons  which  provide  a  potential  well  in  which  the 
quasideuteron  moves.  The  quasideuteron  wavefunction  is  represented 
by  ^(r)  given  by  the  following  expression: 


V  (r) 

k 


(4tt) 


1/2 


sin  (kr+6) 
sin  6 


( a2+k 2)  1/2  v 1/2r 


(1.53) 


where  x  Is  a  function  depending  upon  the  shape  of  the  potential,  v 
the  volume  of  the  nucleus,  and  k  the  wavenumber  corresponding  to  the 
relative  motion  of  the  proton  and  neutron  of  the  quasideuteron. 


k  =  lifj  _  if 2 1 


The  phase  shift  6  is  given  by: 

x  a 

cos  <5  =  -  — 
k 

where  a-1  is  scattering  length.  The  ratio  of  the  cross  section  for 
the  photodisintegration  of  the  quasideuteron  to  the 
photodisintegration  of  the  free  deuteron  is  given  by  the  following 


expression: 


2 7T  (1-otr  ) 
o 


(1.54) 


a(  a^+k  2)v 


where  rQ  is  the  effective  range.  The  cross  section  for  the 
photodisintegration  of  the  nucleus  is  given  by  o: 


(1.55) 


where  the  factor  6.4  is  known  as  the  Levinger  factor. 

This  model  of  the  quasideuteron  was  essentially  a  forerunner 
for  the  Two-Nucleon-Mechanism  (TNM)  or  the  Siegert  term  in  the 
electromagnetic  interaction.  The  TNM  proposed  considered  the 
momentum  transferred  to  the  nucleus  by  the  incoming  particle  as 
being  shared  by  two  nucleons  either  by  the  formation  of  a 
quasideuteron  or  by  A  isobar  production.  These  two  terras  are 
diagramatically  represented  in  fig.  1.4. 

The  TNM  is  less  sensitive  to  the  high  momentum  components  of 
the  wavefunction  since  the  momentum  is  shared  by  the  nucleons  of  the 
p-n  pair.  The  Single-Nucleon-Model  (SNM)  which  could  be  called  the 
direct  knockout  term,  on  the  other  hand,  is  more  sensitive  to  the 
high  momentum  components  of  the  wavefunction,  since  the  momentum 
transferred  is  taken  up  by  a  single  nucleon  in  the  reaction,  say,  P 
+  T  ^He  +  y.  The  one-body  current  term  given  by  eqn.  (1.43) 


stands  for  the  SN  mechanism.  The  differential  cross  section  in  the 
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FIGURE  1.4 

QUAS I- DEUTERON  AND  A-ISOBAR  DIAGRAMS 


TERM  A  :  QUASIDEUTERON 


TERM  B  :  A-ISOBAR 


SNM  is  proportional  to  the  square  of  the  momentum  space  wavefunction 
of  the  bound  proton. 

^  (E  0)  =  C(E  ,0)|  $  (q)|  2  (1.56) 

d  £2  p  p  J ,  L ,  m  j 

where  C  is  the  kinematic  factor  and  <J>  T  is  the  momentum  space 

J,L,mJ 

wavefunction  of  the  nucleon  with  momentum  q  and  specified  by  the 
quantum  numbers  J,  L  and  m^.  The  high  momentum  component  of  the 
nuclear  wavefunction  corresponds  to  the  short  range  correlation  of 
the  nucleons,  dominated  by  the  exchange  of  heavier  mesons  p,  w. 

Ruderman  (Rud-52)  took  up  the  quasideuteron  model  and  got  an 
improved  expression  for  the  cross  section  (see  Chapter  4).  He 
expressed  the  nuclear  wavefunction  in  terras  of  the  wavefunction  for 
the  quasideuteron  and  the  wavefunction  for  the  remaining  nucleons. 

He  evaluated  the  cross  section  for  Pd  -►  t  tt+  in  terras  of  the  cross 
section  for  PP  ■*  d tt4"  on  the  assumption  that  the  neutron  in  the 
triton  does  not  participate  in  the  interaction  and  that  the  incoming 
proton,  after  emitting  a  tt+  forms  a  quasideuteron  with  the  proton  in 
the  deuterium. 

This  model  was  once  more  improved  by  Fearing  (Fea-75a)  (again 
see  Chapter  4)  to  include  distortion  in  the  incoming  channel  and  the 
outgoing  channel.  He  also  used  improved  wavef unctions  for  the 
nuclei.  The  nuclear  cross  section  was  expressed  in  terms  of  the 
cross  section  for  pp->dir+,  a  form  factor  (which  is  the  fourier 
transform  of  the  overlap  of  the  wavef unctions  with  respect  to 
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momentum  transfer),  a  spin  factor  and  a  kinematic  factor.  The  form 
factor  contained  in  it  the  distortion  functions  which  were  derived 
by  integrating  the  distorting  potential  over  the  projectile  path. 

The  distorting  potential,  in  turn,  was  obtained  in  a  Glauber 
formalism  where  NN  total  cross  section  data  and  forward  scattering 
amplitude  data  were  used.  The  distortion  reflected  itself  by  giving 
normalization  to  the  cross  section  whereas  the  form  factor  changed 
the  shape  of  the  angular  distribution  of  the  cross  section.  This  is 
again  because  the  form  factor  is  a  function  of  the  momentum  transfer 
which  is  a  function  of  the  scattering  angle.  The  model  is  known  as 
Distorted  Wave  Impulse  Approximation  (DWIA). 

The  comparison  of  the  experimental  results  of  this  work  with  a 
DWIA  model  calculation  may  be  seen  in  Chapter  4. 

1.5  TESTING  TIME  REVERSAL  INVARIANCE 

This  reaction  provides  a  possible  check  on  Time  Reversal 
Invariance  (TRI)  in  the  electromagnetic  interaction  of  hadrons.  The 
principle  of  detailed  balance  (described  later)  relates  the  cross 
sections  for  the  forward  and  the  inverse  reactions  in  terms  of 
the  momenta  and  spin  multiplicities  of  various  participants  of  the 
reaction.  By  comparing  the  results  of  the  experiment  with  the 
inverse  reaction  data  (Arg-75,  Pic-70,  Are-79,  Sch-83,  Kie-73)  we 
could  determine  the  validity  of  the  time  reversal  principle. 
Described  below  are  some  of  the  aspects  of  TRI  in  the  strong, 
electromagnetic  and  weak  interactions. 

The  discovery  of  parity  violation  in  weak  interactions  in  1957 
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(Hil-58)  triggered  the  search  for  possible  violations  in  strong,  e-m 
and  weak  interactions.  Combined  CP  non-conservation  was  observed  in 
the  decay  pattern  of  the  long-lived  neutral  K  meson  in  1964  by 
Christenson  et  al.  (Chr-64).  The  ratio  of  weak  decay  of  the  into 
two  pions  to  the  sum  of  PC  conserving  decay  rates  was  found  to  be 
close  to  0.2%.  The  overall  conservation  of  CPT  is  expected  to  hold, 
however,  to  a  very  high  precision.  This  follows  from  invariance 
under  relativistic  transformations  in  local  field  theory.  It 
relates  properties  of  particles  and  anti-particles.  The  best 
experimental  verification  is  the  measured  mass  difference  of  K°  and 
K°  where  TCP  is  found  to  hold  at  a  level  of  1  in  lO14.  Thus  we 
conclude  that  CP  invariance  violation  is  associated  T  violation. 
decay  therefore  provides  an  indirect  verification  for  the  TRI  in 
weak  interaction.  This  indirect  T  violation  encouraged  everybody  to 
look  for  a  direct  case  of  T  violation  in  weak,  e-m  and  strong 
interactions.  The  frontiers  of  investigation  are:  checking  the 
principle  of  detailed  balance,  looking  for  T-odd  terras  in  particle 
and  nuclei  decay,  polarisation-analyzing  power  study  for  elastic  and 
inelastic  scattering  and  the  electric  dipole  moment  for  the  neutron. 
Described  below  are  the  details  of  these  aspects. 

It  may  be  noted  that  the  T  operator  is  expressed  as  the  product 
of  a  unitary  operator  and  a  complex  conjugation  operator.  That  is, 

T  is  antiunitary.  This  means  that  there  is  no  conserved  quantity 
associated  with  T  unlike  the  case  of  parity  operator. 

T  invariance  can  be  tested  in  particles  and  nuclei  by  looking 
for  T-odd  terms  which  are  known  as  T  violating  phases.  Let  us 
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consider  the  decay  n  -*•  p  +  e  +  v.  If  a  is  the  neutron  spin  and  p  1 
and  p2  are  the  neutrino  and  electron  momenta,  then  the  terms 
o*(p1xp2)  is  odd  under  time  reversal.  A  non-zero  expectation  value 
of  this  term  is  an  indication  of  T  violation. 

A  number  of  tests  of  T  invariance  in  nuclear  gamma  decays  have 
been  performed  (Gim-79).  Gimlett  et^  al.  looked  at  the  129  keV  gamma 
transition  in  l9jTr.  The  linear  polarization  of  the  129  keV  gamma 
ray  was  measured  directly  with  a  Compton  polarimeter  and  making  use 
of  the  Mossbauer  effect.  The  correlation  terra  (J «kxe) (/*k) ( J «e)  was 
looked  for.  Here  3  is  the  spin  vector  of  the  initial  nuclear  state, 
k  the  photon  momentum  and  e  the  photon  polarization  vector.  This 
term,  as  we  may  see,  changes  sign  under  T  and  so  must  have  zero 
expectation  value  in  the  absence  of  final  state  effects. 

The  absence  of  electric  dipole  moment  for  any  elementary 
particle  is  a  consequence  of  TRI.  In  general,  TRI  prohibits  all  odd 
electric  and  even  magnetic  moments  of  the  elementary  particles.  The 
measurement  of  electric  dipole  moment  (Ram-81)  yields  a  value  of 
|d^|  <  3  x  10"  24  cm,  where  d^  is  the  dipole  moment  per  unit  charge. 
This  may  be  expressed  as  the  product  of  a  parity-violating  coupling 
strength  g^,  a  T-violating  coupling  g^  and  a  hadronic  length 
(Compton  wavelength)  X^.  For  the  neutron  this  product  is 
10~20  gT  cm. 


d 

n 


;Tx 

T  c 


a  5  X  10"  7  X  2  X  10"  14 
=  10"20  gT  cm 
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where  <  4  x  10“ 5,  which  represents  the  most  stringent  limit  for 
T -violation  in  nuclear  physics. 

For  the  elastic  scattering  of  spin  1/2  particles  by  a  target  of 
arbitrary  spin,  polarization  equals  analyzing  power,  P  =  A^.  This 
was  proved  by  Wolfenstein  and  Ashkin  (Wo 1-52).  This  is  a 
consequence  of  the  invariance  of  the  interaction  under  time- 
reversal.  The  P-A  theorem  states  equality  between  the  analyzing 
power  in  the  direct  reaction  and  polarization  in  the  inverse 
reaction  (Sat-58)  for  inelastic  scattering.  This  is  also  a 
consequence  of  TRI.  Note  that  a  non-zero  spin  target  is  necessary, 
otherwise  parity  conservation  alone  dictates  P=A.  The  term  P-A  for 
elastic  scattering  can  be  expressed  in  terms  of  the  product  of  two 
amplitudes,  one  of  which  violates  TRI  while  the  other  is  T- 
invariant.  Measurements  of  P  and  A  for  180  MeV  protons  (Hil-58) 
scattered  from  lithium,  beryllium  and  aluminum  confirmed  P=A  to  a 
few  percent  level.  Comparisons  of  P  and  A  for  p-p  elastic 
scattering  (HUA-6Q)  have  shown  P=A  to  a  few  percent  level. 

Let  us  now  take  a  look  at  one  other  important  aspect  of  TRI, 
namely  the  principle  of  detailed  balance.  The  transition 
probability  from  state  | i>  to  state  |f>  is  related  to  the  S -matrix 
in  the  following  way: 


W±4f  a  |  <f  |  s  |  i> |  2  pf  (1.57) 

where  is  the  density  of  final  states.  The  invariance  of  the 
Hamiltonian  H  under  T  dictates  the  following  equation: 
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ST  =  T  S  r1  -  s*  (1.58) 

where  is  the  time-reversed  S -matrix.  We  get: 

<f|s|i>  =  <iT I s I f T>  (1.59) 

where  i^  and  are  the  time-reversed  states,  that  is,  where  spins 
and  angular  momenta  are  reversed.  The  relation  gives  rise  to  the 
following: 


i*f 


fT  +1T 


(1.60) 


This  equation  yields,  under  the  assumptions  of  parity  conservation, 
absence  of  initial  state  polarization  and  no  final  state 
polarization  detection,  the  following. 


(1-61) 


This  equation  is  known  as  the  principle  of  detailed  balance. 

The  reciprocity  relation  (1.60)  holds  good  for  individual  spin 
states.  Detailed  balance  is  obeyed  when  all  spin  states  are  equally 
populated. 


Let  us  consider  the  reactions  a  +  b  ++  c  +  d.  We  can  convert 


eqn.  (1.61)  to  an  equation  involving  ,  spin  factors  and  centre  of 
mass  momenta,  by  summing  over  final  spin  states,  averaging  over  the 
initial  spin  states,  introducing  the  phase-space  factors  and 
dividing  by  the  incident  particle  flux.  We  get  the  following 
relationship  connecting  the  differential  cross  sections  for  the 
forward  and  the  inverse  reactions. 


(a  +  b  .  c  +  d) 
4?  (c  +  d  -►  a  +  b) 


(2SC  +  l)(2Sd  +  1)  [  Pcd 

T2S  +  1  )(2SU  +  1  J  P  v 

v  a  '* *  v  b  J  ab 


2 

(1.62) 


where  (2S^  +  1)  etc.  are  the  spin  multiplicities,  p^  and  p^  are 
the  CM  momenta.  This  is  also  known  under  the  name  principle  of 
detailed  balance. 

Non-observance  of  the  above  equation  is  an  indication  of  TRI 
violation  whereas  the  observance  of  the  equation  need  not 
necessarily  imply  that  TRI  is  holding  good.  This  is  because  the 
reciprocity  relation  is  guaranteed  by  the  hermiticity  of  the 
Hamiltonian  when  the  transition  probability  is  given  by  Fermi’s 
golden  rule  in  the  first-order  perturbation  theory. 

By  examining  eqn.  (1.62)  it  may  be  seen  that  the  cross  section 
for  the  forward  reaction  is  equal  to  the  cross  section  for  the 
inverse  reaction  times  a  factor  which  is  independent  of  the  angle. 
This  means  that  the  shape  of  the  angular  distribution  for  the 
forward  reaction  is  the  same  as  the  one  for  the  inverse  reaction. 


A  check  of  detailed  balance  was  made  for  the  reaction 


d  +  2^lg  «~»p  +  25Mg  (Wei-68).  It  was  found  to  be  good  within  0.3%. 
Accurate  measurements  of  the  cross  sections  for  the  forward  and 
inverse  reactions  (Dri-79)  have  shown  that  detailed  balance  is 
generally  good  to  0.3%. 

Let  us  now  take  a  look  at  the  detailed  balance  situation  for 
the  electromagnetic  interactions.  Bernstein  et_  al.  (Ber-65) 
suggested  that  the  reactions  of  the  form  y+A^B+Casa 
possible  test  for  TRI.  It  may  be  noted  that  y-N-N  vertex, 
hermiticity  and  current  conservation  obscure  any  T-violation  that 
may  be  present.  Hence  the  above  reactions  must  be  conducted  at 
resonance  energies  so  that  A(1232)  production  replaces  the  y-N-N 
vertex  by  the  y-N-A  vertex. 

Investigation  in  this  regard  has  been  done  for  the  reaction 
p  +  d  «--►  %e  +  y  (Cam-82,  Cam-83)  and  it  is  found  that  there  is  a 
very  good  agreement  in  the  results  which  leads  one  to  believe  that 
TRI  is  holding  good  in  the  e-m  interaction  of  the  hadrons.  It  is 
worthwhile  and  interesting  to  make  another  attempt  in  this  respect 
by  studying  the  reaction  p  +  T  +  ^e  +  y  and  comparing  that  reaction 
to  the  existing  inverse  reaction  data.  This  has  been  one  of  the 
motives  behind  undertaking  to  study  this  reaction. 


CHAPTER  2 


THE  EXPERIMENT 

2.  1  THE  CYCLOTRON 

The  experiment  was  performed  at  the  4BT1  station  of  the  TRIUMF 
Cyclotron  facility  situated  in  Vancouver.  This  is  a  six-sector, 
sector-focussed  negative-ion  cyclotron.  The  beam  extraction  is  done 
by  inserting  a  carbon  foil  to  strip  off  two  electrons  from  the  H“ 
ion.  The  beam  is  continuously  variable  in  energy  from  180  to  520 
MeV,  with  an  intensity  of  100  pA  at  500  MeV.  Polarized  and 
unpolarized  proton  beams  are  available  with  typical  extracted 
intensities  of  300  nA  and  120  pA  respectively.  To  avoid  beam  loss 
by  gas  stripping,  a  vacuum  of  10“ 7  torr  is  maintained.  The 
unpolarized  ion  source  is  of  Ehlers  type,  giving  1-2  inA  of  H“  ions 
at  12  keV.  The  Lamb  Shift  polarized  ion  source  (Don-71)  is  capable 
of  giving  1  pA  of  75-80%  polarized  beam  at  injection.  The 
accelerating  system  consists  of  two  rows  of  quarter-wavelength 
cavities  operating  at  23.055  MHz.  This  provides  a  5  ns  pulse  every 
43  ns. 


2.2  PREPARATORY  CALCULATION 

It  was  important  to  have  an  estimate,  before  the  start  of  the 
experiment,  of  various  things  such  as  the  singles  and  coincidence 
count  rates,  multiple  scattering,  positioning  of  various  detectors 
etc.  It  was  also  decided  to  check  the  working  of  the  detectors  and 
the  electronics  by  measuring  the  cross  section  of  p  +  d  -►  +  y 

(Appendix  D),  which  had  been  measured  earlier  and  therefore  was  well 
known.  These  consistency  checks  would  then  be  a  direct  proof  of  the 
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correct  functioning  of  the  entire  set-up. 

The  target  used  was  TiT,  therefore  the  singles  count  rate  was 
dominantly  determined  by  the  differential  cross-section  for  the 
proton  to  be  scattered  off  Ti.  The  data  was  available  for  the 
cross  sections  -  both  elastic  and  inelastic  -  for  the  Ca  target 
(Wal-66,  Cha-56),  which  is  close  to  that  for  Ti.  The  singles  count 
rate  is  given  by: 


count  =  IN  Aft  2.2.1 

sm  dft 

where  1  is  the  beam  intensity,  N  the  number  of  Ti  atoms  per  cm2  of 
the  target,  ——  the  cross  section  value  at  the  angle  considered  and 

Q  ol 

A  ft  the  solid  angle  subtended  by  the  wire  plane  at  the  target. 

A  typical  value  of  the  singles  count  rate  was  found  to  be  3  x  105/s 

at  an  angle  of  15°  for  300  MeV  protons  with  an  intensity  of  10  nA. 

The  coincidence  count  rate  was  estimated  using  the  inverse 

reaction  data  y  +  4He  -►  P  +  T  (Are-79).  Equating  the  total  CM 

energy  for  the  forward  and  backward  reactions  gave  the  proton  energy 

corresponding  to  the  y  energy.  The  corresponding  cross-section  for 

P  +  T  y  +  ‘♦He  was  obtained  using  the  principle  of  Detailed  Balance 

assuming  Time  Reversal  Invariance.  A  typical  value  of  the 

coincidence  count  rate  for  300  MeV  protons  with  a  10  nA  current  was 

found  to  be  ''■92 /hr  for  =  60°  and  ~1  /hr  for  =  150°. 

Y  Y 

The  angular  spreading  of  ^He  was  calculated  from  the  knowledge 
of  beam  spreading  and  multiple  scattering.  This  calculation  was 


important  in  determining  the  distances  of  the  front  detectors  from 
the  target  by  making  sure  that  the  4He  spreading  did  not  exceed  the 
boundaries  of  the  Cherenkov  counter  image  on  the  wire  plane. 

2.3  EXPERIMENTAL  SET-UP 

The  experiment  was  carried  out  in  a  two-arm  coincidence  mode. 
The  recoil  ^He ’ s  were  detected  in  coincidence  with  the  corresponding 
Y  on  the  opposite  side.  The  detectors  were  arranged  on  either  side 
of  the  beam  so  that  we  had  two-fold  coincidence  measurement 
(Fig.  2.1).  By  convention,  the  side  was  determined  looking 
downstream  from  the  target. 

The  cross  section  and  the  analyzing  power  were  studied  at  three 

different  proton  energies;  227  MeV,  300  MeV  and  375  MeV.  The 

measurements  were  made  at  0^,M  -  90°  for  227  MeV,  =  54°,  65°, 

78°,  90°,  99°,  118°  for  300  MeV,  and  0^M  =  60°,  91°  for  375  MeV. 

These  provided  enough  data  to  compare  with  the  inverse  reaction, 

thereby  enabling  us  to  draw  conclusions  on  the  validity  of  Time 

Reversal  Invariance.  The  overall  integrated  charge  at  300  MeV  for 

0^M  -  65°,  90°,  118°  was  roughly  415  pc  for  the  target  run  and 

230  pc  for  the  background  run.  A  total  of  2.2  M  events  were 

recorded,  17,000  due  to  the  pulser.  The  events  due  to  (p,y)  were 

found  to  be  1502  at  0^  =  65°,  441  at  90°,  287  at  118°,  counts 

included  for  both  sides.  The  corresponding  numbers  for  the 

background  run  were  0.1  M,  the  pulser  contributing  20,000.  For  the 
CM 

same  energy  at  0^  =  54°,  78°,  99°  the  integrated  charge  was  875  pc 
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(1.1  M  events,  0. 2  M  pulser  events)  for  the  target  run  and  446  jjc 
(0.6  M  events,  60,000  pulser)  for  the  background  run.  The  (p,y) 
events  were  3601  at  54°,  1770  at  78°  and  860  at  99°.  A  total  of 
350  yc  was  collected  for  227  MeV,  corresponding  to  0.5  M  events  and 
22,000  pulser.  The  background  run  had  97  yc  corresponding  to  0.1  M 
events  and  2400  pulser  events.  560  events  corresponded  to  (p,y)  at 
90°CM.  The  integrated  charges  at  375  MeV  were  206  yc  (0.6  M  events, 
16,000  pulser)  and  63  yc  (0.1  M  events,  5000  pulser)  for  the  target 
run  and  the  background  run  respectively.  (p,y)  events  were  816  at 
60°  and  285  at  91°. 

It  may  be  noted  that  the  finite  angular  acceptance  provided  by 
the  recoil  particle  detector  made  it  possible  to  detect  y’s  at  three 
different  angles  simultaneously.  A  typical  value  of  the  angular 
acceptance  of  the  front  detector  was  ±4°.  At  300  MeV,  an  angular 
coverage  of  10.8°  to  18.8°  enabled  measurements  of  y  at  angles  (lab 
frame)  50°  to  60°,  73°  to  83°  and  102°  to  112°,  which  corresponded 
to  254-262  MeV,  234-243  MeV  and  213-219  MeV  y's.  The  Cherenkov 
counters  themselves  had  an  angular  acceptance  of  ±5°. 

The  entire  experiment  involving  the  above  mentioned 
measurements  along  with  one  run  for  p  +  d  %e  +  y  took  up  74  runs 
spanned  over  a  period  of  10  days. 

2.4  THE  SCATTERING  CHAMBER 

A  rectangular  evacuated  chamber  of  dimensions  33"  x  23  1/2"  x 
19"  was  used  to  house  the  target.  The  chamber  had  two  front  windows 
separated  by  3  1/4".  Each  window  was  of  dimensions  9  1/4"  x  3". 

The  side  windows,  one  on  either  side  of  the  chamber,  were  of 
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dimensions  32"  x  5".  The  back  windows,  separated  by  3  1/4",  were  of 
dimensions  9  1/4"  x  3".  Kapton  of  thickness  127  was  used  for  the 
windows.  The  dimensions  of  the  front  windows  provided  an  angular 
sweep  of  approximately  5.5°  to  33°  on  either  side  for  the  recoil 
particle.  The  side  window  dimension  allowed  an  angular  sweep  of  36° 
to  144°  for  photons  on  either  side.  The  back  windows  allowed 
measurements  up  to  174°. 

2.5  THE  TARGET 

Tritiated  titanium  foil  with  a  copper  backing  was  used  as  the 
target.  Background  counts  were  obtained  by  using  an  identical 
nontritiated  foil.  The  target,  supplied  by  ORNL,  was  made  by 
absorbing  tritium  into  the  crystalline  structure  of  titanium 
evaporated  onto  a  copper  backing.  When  in  use  the  target  was  placed 
in  an  aluminum  enclosure  maintained  at  a  vacuum  of  100  \i  torr.  2"  x 
2"  tantalum  squares  sandwiched  the  titanium  foil.  A  0.78"  diameter 
hole  at  the  center  of  the  tantalum  square  exposed  the  tritiated 
region.  This  was  held  at  the  center  of  the  aluminum  enclosure  by 
means  of  foam  rubber  spacers  at  the  corners  compressed  between  two 
window  frames.  A  Kapton  window  of  1  mil  thickness  was  glued  onto 
the  inside  surface  of  each  window  frame.  The  frame  was  sealed  onto 
the  body  of  the  enclosure  by  means  of  "0"  rings. 

To  determine  the  tritium  content,  the  cross  section  for  (P,T) 
elastic  scattering  was  measured  using  the  MRS  (Appendix  A).  Taking 
into  account  the  decay  loss,  the  thickness  was  deduced  to  be  1.28 
mg/cm2  at  the  time  the  experiment  was  done.  The  other  constituents 
of  the  target  were  Ti  (19.4  mg/cm2)  and  Cu  (5  mg/cm2). 
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The  targets  mounted  on  the  target  ladder  were  tritiated 
titanium  foil,  nontritiated  titanium  foil,  a  ZnS  screen  and  a  CD^ 
target.  The  ZnS  helped  locate  and  center  the  beam. 

2. 6  THE  HELIUM  BAG 

The  space  between  the  scattering  chamber  front  window  and  the 
wire  chamber  plane  was  filled  by  a  plastic  bag  holding  helium  at 
atmospheric  pressure.  This  helped  reduce  the  multiple  scattering  of 
the  recoil  particles  on  the  way  to  the  detector. 

2.1  THE  6-RAY  MAGNET 

A  description  of  6-rays  is  given  in  Appendix  B.  It  was 
important  to  prevent  6-rays  from  getting  into  the  front  detectors. 

A  magnet  of  5  kG  field  strength  was  used  to  deflect  away  the  6-rays. 
The  field  was  vertical  so  that  6-rays  were  deflected  towards  the 
right.  The  presence  of  the  magnet  was  taken  into  account  and  the 
angles  of  the  wire  chambers  were  reset  so  that  they  would  again  be 
centered  on  the  recoil  particles  scattered  at  the  kinematically 
correct  angles.  An  offset  of  approximately  1°  was  required  for  each 
wire  plane. 

2.8  THE  POLARIMETER 

The  beam  current  and  polarization  were  measured  by  an  in-beam 
polarimeter  situated  at  a  distance  of  250  cm  upstream  from  the 
target.  A  detailed  description  appears  elsewhere  (Gre-79). 

Fig.  2.2  shows  the  polarimeter  set  up.  The  target  used  was  a  CH 2 
foil  of  thickness  5.57  mg/cm2.  The  polarimeter  monitored  the  pp 
elastic  scattering  and  the  incident  proton  flux  was  elucidated  from 
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the  known  cross  section  at  O^ab  =  ^°*  Beam  polarization  was 
determined  from  the  p-p  analyzing  power,  which  is  also  well  known. 
The  value  of  the  analyzing  power  at  Q^ab  =  is  at  300  MeV 

incident  proton  energy. 

Plastic  scintillator  detectors  LI,  L2  and  R1 ,  R2  were  placed  at 
a  lab  angle  of  17°.  The  corresponding  recoil  protons  were  detected 
by  the  scintillators  R3  and  L3.  Because  of  the  low  energy  of  the 
recoil  protons,  backward  detectors  were  situated  inside  the 
polarimeter  scattering  chamber.  LI  and  R1  determined  the  solid 
angle.  These  forward  angle  defining  detectors  were  rotated  such 
that  normals  drawn  from  their  centers  coincided  at  a  point  on  the 
beam  downstream  from  the  polarimeter  target.  This  positioning 
helped  minimize  differences  in  the  left  and  right  solid  angles  due 
to  non-centering  of  the  beam  on  the  target.  The  effect  was  to 
increase  (decrease)  the  solid  angle  at  larger  (smaller)  scattering 
angles  so  that  the  count  rate  remained  the  same. 

Fig.  2.3  shows  the  coincidence  circuit  associated  with  the 
polarimeter.  The  event  was  determined  by  LI *L2 «R3  (left)  and  by 
R1 *R2 *L3  (right).  Background  from  carbon  was  independently 
determined  using  a  carbon  target.  L1«L2*R3  gave  the  contribution 
due  to  reals  and  accidentals  for  the  left  side.  The  accidental 
contribution  alone  was  determined  by  a  coincidence  rate  of  R3  with 
L1»L2  delayed  by  43  ns  (the  separation  between  consecutive  beam 
bursts) . 

No.  of  incident  protons  =  (  — — ^-c-c^  -r — — — -  )  (2.8.1) 
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where  P£,  PJiacc  refer  to  the  total  counts  and  accidental  counts 
respectively  for  the  left  side:  F  is  the  polarimeter  factor  which 

stands  for  counts  per  nanocoulomb. 

,  P  -P„  -P  +P 

ti  1  /  i  cc  r  race  x  „  ON 

Beam  polarization  =  —  (  p—p - +p  _.p -  )  (2.8.2) 

i  Jiacc  r  race 

where  A  is  the  polarimeter  analyzing  power.  The  contribution  to  the 
analyzing  power  from  C(p,2p)  was  accounted  for  by  using  the  carbon 
target. 

2.9  THE  DETECTORS 

The  a  particles  were  detected  in  an  E-AE  counter  telescope. 
Looking  downstream,  this  counter  assembly  followed  a  multiwire 
proportional  chamber.  The  y  rays  were  detected  in  lead  glass 
Cherenkov  counters.  These  were  preceded  by  a  veto  counter  to 
discard  events  arising  from  charged  particles  (pions  or  electrons) 
getting  into  the  Cherenkov  counter.  A  detailed  description  of 
various  detectors  follows. 

2.9.1  Charged  Particle  Counter  Telescope 

The  AE  detector  was  a  plastic  scintillator  with  dimensions 
1/25"  x  5"  x  5”.  This  was  coupled  to  an  RCA  8575  photomultiplier 
tube  through  a  light  pipe.  The  E  detector  was  a  plastic 
scintillator  with  dimensions  1/2”  x  5"  x  5",  again  coupled  to  an  RCA 
8575  photomultiplier  tube  through  a  light  pipe.  The  scintillator 
material  was  NE110  plastic.  The  plastic  was  chosen  to  give 


reasonable  energy  resolution  and  to  accommodate  the  high  count  rate 
expected  in  these  detectors  (see  Sec,  2.2),  NE110  has  good  light 

transmission  and  a  decay  time  of  3.3  ns. 

2.9.2  Multiwire  Proportional  Chamber 

Multiwire  proportional  chambers  were  used  to  obtain  positional 
information  on  the  recoil  particles.  Each  chamber  was  5"  x  5" 
having  an  x  plane  and  a  y  plane  of  64  wires.  Details  of  these  MWPC 
are  described  elsewhere  (Cai-75) .  The  chambers  had  a  -5  kV  high 
voltage  plane,  an  x  sense  plane,  another  high  voltage  plane,  a  y 
sense  plane  and  yet  another  high  voltage  plane.  The  high  voltage 
wires  were  102  ym  in  diameter  and  were  spaced  1  mm  apart.  The  sense 
plane  wires  had  separations  of  2  mm  and  were  16  yn  in  diameter.  The 
sense  plane  wires  were  of  gold-plated  tungsten  whereas  the  HV  plane 
wires  were  of  beryllium  and  copper.  The  chambers  had  25  \m  thick 
kapton  windows.  A  magic  gas  mixture  was  circulated  through  the 
chamber.  The  constitution  of  the  magic  gas  mixture  was  75%  Ar 
(which  was  bubbled  through  methylal  CH  2  (OCH  3) 2  at  0°C)  ),  24.5% 
isobutane  C^H^,  and  0.5%  freon  -  13B1  (CBrF3).  For  details  refer 
to  (Cha-72).  Argon  provides  the  amplification,  isobutane  absorbs  uv 
radiation  (thus  extending  the  proportional  region  allowing  higher 
amplification),  freon  captures  electrons  (thus  providing  heavy 
negative  ions  near  the  positive  ions  and  avoiding  space-charge 
effects).  It  also  avoids  avalances.  Methylal  eliminates  the 
degradation  of  MWPC 1 s  with  integrated  particle  fluxes.  The  wire 
chamber  readout  system  was  of  the  amplifier  per  wire  type  (Cai-75). 


2.9.3  The  Cherenkov  VETO  Detector 


A  plastic  scintillator  (NE110),  6"  x  6”  x  1/8",  was  used  as  a 
veto  counter  before  each  Cherenkov.  This  was  coupled  to  an  RCA  8573 
photomultiplier  tube  through  a  light  pipe.  Charged  particles 
entering  the  Cherenkov  counter  which  could  cause  false  events  were 
eliminated  using  signals  from  the  veto’s. 

2.9.4  The  Cherenkov  Counter 

Lead  glass  Cherenkov  counters  7"  long  x  5"  in  diameter  were 
used  to  detect  y’s  (see  fig.  2.4).  These  were  coupled  to  an 
RCA  8854  photomultiplier  tube  with  RTV  602  silicon  cement.  The 
composition  of  lead  glass  by  weight  was:  Pb,  51%;  0,  26%;  Si,  18%, 
K,  3%;  and  Na,  1.5%.  The  y's  in  the  lead  glass  undergo  two  dominant 
mechanisms:  pair  production  and  compton  scattering.  At  higher 

energies  (a  typical  value  in  the  present  situation,  for  the  y  energy 
is  270  MeV  for  300  Mev  proton  energy)  pair  production  dominates. 

The  resulting  electromagnetic  cascade  shower  yielded  enough  light  to 
be  seen  by  the  photomultiplier  tube. 

Since  the  plastic  scintillator  had  enough  resolution  (6%)  to 
separate  the  a’s  corresponding  to  (p,tt°)  and  (p,y)  reactions,  it  was 
not  necessary  to  have  a  high  resolution  y  detector.  The  lead  glass 
detectors  gave  y  detection  with  high  efficiency  (~90%). 

The  efficiency  of  the  Cherenkov  counters  was  estimated  using  a 
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Monte-Carlo  technique  based  on  the  one  described  by  Sober  (Sob-73). 
Electron  calibration  data  (Appendix  C)  was  used  in  this  calculation 
to  determine  Cherenkov  counter  efficiency  to  detect  y’s. 

Fast  signals  from  the  E,  AE  and  C  detectors  were  used  to  form 
a  trigger.  The  Multiwire  Proportional  Chamber  (MWPC)  provided  the 
positional  information  for  ^He  without  participating  in  the  event 
defining  logic. 

2.10  THE  ELECTRONICS 

The  fast  electronics  used  to  form  the  trigger  (figs.  2.5,  2.6, 
2.7)  were  located  in  the  experimental  hall.  This  part  of  the 
electronics  was  used  in  the  event  formation  and  multiwire 
information  processing.  Fig.  2.5  gives  the  left -hand -side  fast 
electronics.  The  photomultiplier  signal  pulse  was  split  into  two 
portions.  One  portion  was  fed  into  the  Analogue  to  Digital 
Converter  (ADC).  The  ADC  unit  used  was  LRS  CAMAC  model  2249.  This 
is  a  high  resolution  (0.1%),  multi-input  (12-fold)  ADC  of  charge 
integrating  type.  Input  sensitivity  is  0.25  pc/count.  This  has  a 
dynamic  range  of  1024  channels.  The  other  portion  of  the  signal  was 
fed  into  a  discriminator.  A  constant  fraction  discriminator  CFD  934 
(EG  &  G  ORTEC  934)  was  attached  to  E,  AE  and  C  detectors  whereas  the 
LRS  NIM  Model  621  was  attached  to  the  veto.  CFD  934  is  a  quad 
discriminator  with  a  double-pulse  resolution  of  <  10  ns  and 
threshold  varying  from  -30  to  -1000  mV.  This  module  splits  the 
input  signal  into  two,  delays  one  portion,  and  subtracts  that  from  a 
constant  fraction  of  the  undelayed  signal.  The  resulting  signal  has 
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a  baseline  crossover  independent  of  the  input  signal  amplitude. 

This  zero-crossing  line  is  used  to  provide  a  precisely  timed  logic 
pulse.  The  output  pulse  is  of  width  6  ns  to  150  ns  with  an 
amplitude  of  -800  mV.  LRS  621  is  a  quad  discriminator  with  a 
threshold  of  -30  mV  to  -1  V  and  an  output  of  -700  mV  to  -800  mV  with 
a  minimum  width  of  5  ns.  The  maximum  pulse  rate  that  can  be  handled 
is  110  MHz. 

One  of  the  outputs  from  the  discriminator  is  fed  into  a 
coincidence  unit  (LRS  NIM  Model  365).  The  other  output  goes  to  a 
Time  to  Digital  Converter  (TDC)  (LeCroy  CAMAC  Model  2228A).  The 
third  output  goes  to  a  scaler  (LRS  CAMAC  Model  84).  Coincidence 
unit  365  has  a  double  pulse  resolution  of  5  ns.  The  LeCroy  2228A  is 
an  octal  TDC  with  an  11  bit  resolution  out  of  2048  channels.  This 
measures  the  time  from  the  leading  edge  of  the  common  start  pulse  to 
the  leading  edge  of  the  individual  stop  pulse.  This  had  100,  200, 
and  500  ns  full-scale  time  ranges  and  provides  50,  100  and  250 
ps/count  resolution  respectively.  Each  channel  disregards  stop 
signals  received  before  start  signals  and  accepts  only  one  stop  for 
one  start.  The  LRS  Model  84  scaler  provides  24  bits  of  data  storage 
which  may  be  stored  at  a  rate  of  100  MHz.  Maximum  capacity  is 
approximately  17  M  counts. 

The  outputs  from  the  three  AND  gates  were  fed  into  an  OR  gate 
and  its  output  defined  EVENT  FAST  for  that  side.  The  outputs  from 
the  two  sides  were  fed  into  an  OR  gate  and  its  output  was  fed  into 
the  LATCH  (LeCroy  NIM  Model  364)  preceding  the  MWPC  controller. 

LATCH  caused  the  wire  chamber  information  to  be  held  until  the  slow 


electronics  determined  whether  the  event  was  good.  The  MASTER  GATE 
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(MG)  signal  (fig.  2.9),  among  other  things,  caused  the  wire  chamber 
information  to  be  read  into  the  computer  memory.  The  data  was  read 
in  via  the  CAMAC  INPUT  REGISTER  (CAMAC  Model  IR-1).  Once  the 
information  was  read  the  controller  was  reset  with  a  CLEAR  signal, 
initiated  by  an  ’EVENT -DONE'  signal  generated  by  the  computer 
through  OUTPUT  REGISTER  5  (Fig.  2.9).  This  CLEAR  signal  also 
unlocked  the  LATCH  imposed  by  EVENT  FAST,  thereby  readying  the  LATCH 
for  another  event. 

Fig.  2.8  shows  one  section  of  the  SLOW  ELECTRONICS.  The  LEFT 
EVENT  FAST  was  set  in  coincidence  with  the  signal  from  LEFT  E. 
Similarly  the  RIGHT  EVENT  FAST  was  set  in  coincidence  with  the 
signal  from  RIGHT  E.  The  outputs  were  fed  into  an  OR  gate  (LRS 
Model  428)  and  subsequently  fed  into  the  MG.  Note  that  LEFT  EVENT 
FAST  as  well  as  RIGHT  EVENT  FAST  were  set  into  coincidence  with  the 
PULSER  signal  generated  by  the  LRS  Model  8010  PULSE  GENERATOR.  This 
was  so  that  the  MG  would  accept  the  PULSER  events  also,  thus  making 
it  available  for  dead-time  determination.  One  output  of  the  pulser 
fan-out  was  fed  into  the  digital  coincidence  register  (DCR)  (CAMAC 
C212  Unit)  which  generated  a  24-bit  word,  each  bit  representing  the 
status  of  an  input. 

Fig.  2.9  gives  the  final  phase  of  the  electronics.  The  PRE¬ 
EVENT  started  the  LEFT  TDC  and  RIGHT  TDC.  This  also  gated  the  FRONT 
ADC.  The  output  signal  of  a  Cherenkov  counter  in  coincidence  with 
the  PRE-EVENT  gated  the  BACK  ADC.  The  PRE-EVENT  also  served  as 
input  to  the  MG.  Output  from  the  MG  initiated  storage  of  the  event 
structure  into  the  memory.  This  also  strobed  the  DCR.  Once  the 
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event  was  written,  the  computer  would  issue  a  CLEAR  signal  through 
OUTPUT  REGISTER  5.  This  reset  the  CAMAC  units  ADC,  TDC,  DCR  and 
also  cleared  the  LATCH  attached  to  the  MWPC  controller.  While  the 
computer  was  busy,  the  MG  was  latched  so  that  no  other  signal 
could  form  a  MG  during  that  time. 

The  LATCH  is  associated  with  the  MWPC  controller  and  the  MG 
accounted  for  the  electronic  dead-time  of  the  system.  The  rest  of 
the  dead-time  was  taken  up  by  the  computer  processing  the  events. 

The  dead-time  of  the  system  was  inferred  from  the  pulser,  presented 
to  the  system  (known  from  the  pulser  scaler),  and  the  pulser 
accepted  by  the  system  (known  from  DCR  information). 

2.11  THE  CAMAC 

The  CAMAC,  standing  for  Computer  Automated  Measurement  and 
Control,  is  an  instrumentation  system  developed  to  accommodate 
digital  circuits.  A  multiwire  printed  circuit  board  is  mounted  at 
the  rear  of  the  power  crate.  This  provides  interconnections  between 
the  modules  and  is  known  as  CAMAC  DATAWAY.  CAMAC  DATAWAY  allows 
bidirectional  communication  between  modules  and  computer  or  between 
modules.  The  crate  used  in  the  experiment  had  25  stations  (slots) 
to  accommodate  the  CAMAC  modules.  Each  slot  had  an  86  pair  card- 
edge  connector  by  which  the  modules  gained  access  to  the  wiring  of 
the  crate  (dataway).  Control  of  modules  within  the  crate  was 
provided  by  the  standard  type-A  crate  controller  which  accepts 
external  commands  from  the  branch  highway  (path  between  crates  and 
computer).  The  rightmost  slots,  i.e.  the  24th  and  25th  slots  were 


occupied  by  a  3900  type-A  crate  controller.  The  rest  of  the  slots 
were  occupied  by  CAMAC  units  such  as  the  ADC,  TDC,  DCR,  etc. 

Figs.  2.10,  2.11  and  2.12  show  the  allocation  of  various  slots 
for  the  CAMAC  modules  used  in  the  experiment.  These  modules  were 


spread  over  three  crates. 


CHAPTER  3:  THE  DATA  ANALYSIS 


3.1  THE  DATA  ACQUISITION 

The  data  was  collected  on  magnetic  tape,  using  a  Honeywell  316 
computer,  on  an  event  by  event  basis.  The  data  contained  the  good 
events  and  the  pulser  events.  The  acquisition  program,  'HONDA',  is 
described  in  detail  in  Hut-81.  The  program  'HONDA'  controlled  the 
3-crate  CAMAC  branch  by  means  of  a  Daresbury  DNPL-EC321A  interface. 
In  addition  to  writing  the  data  onto  the  tape,  the  program  also 
stored  selected  events  in  the  computer  memory,  thus  facilitating  the 
display  of  various  spectra  on  the  printer  or  storage  scope.  Events 
to  be  written  on  tape  could  be  subjected  to  certain  criteria  and 
software  cuts  determined  from  the  displayed  spectra. 

3.2  THE  DATA  STRUCTURE 

The  data  was  written  onto  the  tape  in  blocks  of  length  516 
words,  each  word  made  up  of  16  bits  (Rog-77,  Hut-81).  The  first 
4  words  contained  information  used  in  rereading  of  the  tapes  and  the 
remaining  512  words  contained  the  data.  Each  run  recorded  had  a 
beginning  block  and  an  end  block,  giving  information  on  the  various 
parameters  scaled  during  the  run.  These  two  blocks  sandwiched  a 
variable  number  of  blocks  each  containing  a  variable  number  of 
records,  a  record  representing  the  various  parameters  such  as  ADC, 
TDC,  DCR  and  wire  chamber  co-ordinates,  pertaining  to  a  particular 
event.  Table  3.1  lists  the  scalers  relevant  to  a  run.  Table  3.2 
lists  parameters  representing  an  event.  Since  the  number  of  wires 


73 


74 


TABLE  3.1 
SCALERS 


No. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
lb 

17 

18 

19 

20 
21 
22 

23 

24 

25 
2b 

27 

28 

29 

30 


Scaler 


MASTER  GATE 
CR1  +  CR2  +  CR3 
CL1  +  CL2  +  CL3 
LF  •  RB 


RF 

EL 

ER 

EL 

ER 


LB 


LF 

RF 


RB 

LB 


PULSER 

SECONDARY  EMISSION  MONITOR 
POLARIMETER  LEFT 
POLAR IMETER  LEFT  ACCIDENTALS 
POLARIMETER  RIGHT 
POLARIMETER  RIGHT  ACCIDENTALS 
AEL 


AER 
CL  l 
CL2 
CL3 
CR1 
CR2 
CR3 
VL1 
VL2 
VL3 
VR1 
VR2 
VR3 

THE  ELAPSED  TIME 


LEGEND 

CR1 : 

Cherenkov  Right  1 

CL1: 

Cherenkov  Left  1 

LF  »RB : 

Left  Front  •  Right 

Back 

RF  *LB: 

Right  Front  •  Left 

Back 

EL: 

Left  E 

AEL: 

Left  AE 

ER: 

Right  E 

AER: 

Right  AE 

VL1 : 

Veto  Left  1 

VR1 : 

Veto  Right  1 

TABLE  3.2 


THE  EVENT  STRUCTURE 


Word 

No. 

Parameter 

0 

EVENT  LENGTH 

1 

FLAG 

2 

LAM  PATTERN 

3-4 

DCR 

5 

EL  — | 

6 

CR1 

7 

CR2 

8 

CR3 

TDC  LEFT 

9 

VR1 

10 

VR2 

11 

VR3  — 1 

12 

ER  — | 

13 

CL1 

14 

CL2 

15 

CL3 

TDC  RIGHT 

16 

VL1 

17 

VL2 

18 

VL3  - 

19 

ML  - 

20 

EL 

ADC  FRONT 

21 

MR 

22 

ER  - 

23 

CL1  — j 

24 

CL2 

25 

CL3 

26 

CR1 

27 

CR2 

28 

CR3 

ADC  BACK 

29 

VL1 

30 

VL2 

31 

VL3 

32 

VR1 

33 

VR2 

34 

VR3  — 

35 

• 

— 

MWPC  DATA 

• 

• 
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hit  is  quite  arbitrary,  the  event  length  is  a  variable.  The  ’FLAG' 
in  the  event  structure  is  to  indicate  whether  it  is  the  last  event 
or  not.  The  fLAM  PATTERN'  is  the  CAMAC  'look-at-me'  signal  issued 
by  each  CAMAC  module  when  it  needed  the  attention  of  the  computer. 
'DCR'  stands  for  Digital  Coincidence  Register.  This  is  a  24-bit 
word  of  which  only  8  bits  were  used  in  the  experiment.  Each  bit 
represents  a  characteristic  of  the  event.  Table  3.3  gives  the  DCR 
bit  pattern. 

3.3  THE  DATA  RE-ORGANIZATION 

The  total  data  collected  during  the  experiment  was  recorded  on 
63  tapes  each  of  2400',  at  800  bpi.  This  was  quite  an  intimidating 
size  and  to  facilitate  the  analysis  it  was  decided  to  reorganize  the 
data  to  make  it  more  compact.  The  criteria  used  to  achieve  this  are 
described  below. 

The  information  on  the  BUSY  state  of  the  beam,  available  from 
the  ion  source  spin  indicator,  was  used  to  reject  events  for  which 
the  beam  spin  was  not  clearly  defined.  These  events,  nevertheless, 
were  recorded  on  the  tape.  They  were  discarded  when  the  data  was 
rewritten. 

The  MWPC  data  was  of  variable  length,  comprised  of  each  wire 
with  a  detectable  charge  collected  for  both  the  planes.  Only  the 
centroid  of  the  wires  hit  was  important  in  the  recoil  angle 
calculation  and  so,  the  words  standing  for  the  addresses  of  the 
wires  hit  were  reduced  to  a  single  word  representing  the  centroid  of 
the  wires  hit.  Another  word  indicating  whether  there  was  a  'miss' 


TABLE  3.3 
DCR  BIT-PATTERN 


DCR 

Bit 

Parameter 

0 

LF  •  RB  •  LE  (LEFT  EVENT) 

1 

RF  •  LB  •  RE  (RIGHT  EVENT) 

2 

LF  •  RB  (LEFT  EVENT;  E  not  in  trigger) 

3 

RF  •  LB  (RIGHT  EVENT;  E  not  in  trigger) 

6 

PULSER 

7 

BEAM  SPIN  UP 

8 

BEAM  SPIN  DOWN 

9 

BEAM  SPIN  OFF 

or  a  'multiple’  was  also  included.  A  'miss'  occurred  when  no  wire 
was  hit.  A  'multiple*  occurred  when  the  separation  between  the 
consecutively  hit  wires  was  greater  than  or  equal  to  3.  Events 
characterized  by  a  'miss*  or  'multiple*,  denoted  by  'MM',  were 
discarded  when  scatter  plots  or  histograms  were  generated.  They 
thus  represented  an  inefficiency  in  the  system  and  an  appropriate 
correction  was  applied.  To  sum  up,  wire  chamber  data  took  up  4 
words,  2  for  each  plane  on  the  rewritten  tapes. 

The  large  flux  of  elastically  scattered  protons  which 
infiltrated  the  recoil  particle  detector  assembly  caused  spurious 
events  which  took  up  a  good  portion  of  the  tape.  A  hardware  attempt 
was  made  to  cut  these  down  to  a  manageable  level  by  boosting  the 
threshold  of  a  Constant  Fraction  Discriminator  attached  to  the  AE 
detector.  Some  of  the  remaining  protons  were  tracked  down  through  a 
software  cut  involving  the  sum  of  the  pulse  heights  in  the  AE 
plastic  and  E  plastic.  These  protons,  however,  left  the  leading 
four  words  of  the  event  structure  on  the  tape.  These  events,  called 
'BAD'  on-line,  were  discarded  while  the  data  was  rewritten. 

As  may  be  seen  from  table  3.2,  the  event  structure  contained 
possible  information  from  three  different  gamma  detectors  which  was 
the  general  case.  A  particular  event  would  however  contain  only 
data  from  one  detector  at  a  time.  This  meant  it  was  possible  to 
discard  the  words  for  the  two  non-struck  detectors  which  were 
present  in  the  original  data  sructure.  This  could  be  achieved  if  it 
was  known  which  gamma  detector  gave  rise  to  the  coincidence.  A  cut 
on  the  Cherenkov  counter  TDC  was  imposed  to  make  such  an 


identification.  By  examining  the  TDC  spectrum,  it  was  found  that  an 
event  falling  between  TDC  channels  15  and  645,  which  correspond 
roughly  to  135  ns,  would  represent  a  signal  in  that  detector.  This 
identification  was  called  time -of -flight  test  (TOF  test).  We  thus 
ended  up  with  ten  words  per  event,  as  shown  in  table  3.4.  TOF  TEST 
WORD  -  1,  2,  3,  6,  11,  12,  13  stood  for  CR1,  CR2,  CR3,  pulser,  CL1, 
CL2  and  CL3  events  respectively.  Y  MM  *  1  stood  for  the  fact  that 
the  y  plane  had  a  ’multiple'  or  'miss'.  X  MM  =  1  was  the 
corresponding  word  for  the  x  plane.  Identification  of  the  event 
this  way  had  an  additional  advantage,  since  the  most  forward 

Cherenkov  counters  situated  at  6  *  54°  and  65°  in  two  different 

c .  m. 

configurations,  corresponded  to  57  MeV  and  65  MeV  a's  (for  Ep  » 

300  MeV)  which  stopped  either  in  the  AE  plastic  or  in  the  wrappings 
between  AE  plastic  and  E  plastic  and  hence  the  corresponding  event 
could  not  be  categorized  LEFT  or  RIGHT  since  the  LEFT  or  RIGHT  DCR 
bit  was  defined  with  E  in  coincidence.  Note  also  that  the  original 
DCR  word  was  retained  in  order  to  identify  beam  polarization  as  well 
as  to  retain  information  on  LEFT  or  RIGHT,  if  necessary. 

The  data  was  summarized  into  eleven  2400',  1600  bpi  tapes, 
enabling  an  accelerated  pace  in  further  data  analysis. 

3.4  VARIOUS  STAGES  OF  DATA  ANALYSIS 

3.4.1  Recoil  Angle  Calculation 

The  recoil  angle  was  calculated  from  the  x-plane  centroid  x  and 
the  y-plane  centroid  y.  Fig.  3.1  shows  the  wire  plane  ABMD  situated 
at  a  distance  *d '  from  the  target  T.  The  line  TO  represents  the 
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TABLE  3.4 


REWRITTEN  EVENT  STRUCTURE 


Word 

No. 

Parameter  j 

1 

ORIGINAL  OCR  WORD 

2 

TOF  TEST  WORD 

3 

CHERENKOV  COUNTER  TDC 

1  4 

CHERENKOV  COUNTER  ADC 

5 

AE  PULSE  HEIGHT 

6 

E  PULSE  HEIGHT 

7 

Y  MM  — 

8 

Y  CENTROID 

MWPC  Information 

9 

X  MM 

10 

X  CENTROID  ~ 

beam  line.  C  is  the  center  of  the  wire  plane  and  OC  is  the  normal 


drawn  to  line  TC,  along  the  wire  plane  <OTC  *  0,  the  angle  at  which 
the  wire  plane  is  placed.  The  angle  made  to  the  beam  line  by  the 
line  joining  the  target  T  and  any  point  on  the  wire  plane,  P  is 
denoted  by  6.  PX  is  the  normal  erected  from  P  to  OC. 


3.4. 1.1 


where  OT  *  d  sec  0 


where  and  y^  stand  for  wire  numbers  at  point  c,  and  f  the 
interwire  distance  (2  mm). 

3.4.2  Combining  E  and  AE 

The  problem  of  the  a’s  corresponding  to  the  most  forward 
Cherenkov  counter  stopping  in  the  AE  scintillator  necessitated 
combining  of  the  pulse  height  AP  in  the  AE  plastic  and  the  pulse 
height  P  in  the  E  plastic  with  appropriate  weights  to  obtain  the 
total  energy  ETOT  of  the  a.  The  relative  gains  of  the  E  detector 
and  the  AE  detector  were  necessary  to  add  P  and  AP.  This  was 
obtained  by  looking  at  the  AP  v/s  P  scatterplot  of  the  passing 
protons  (fig.  3.2)  and  also  by  knowing  the  thicknesses  of  the  E  and 
the  AE  detectors.  Passing  protons  are  represented  by  the  straight 
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E  PLASTIC  PULSE  HEIGHT 


line  portion.  The  energy  loss  for  these,  to  a  good  approximation, 
is  directly  proportional  to  the  thickness  of  the  plastic.  Let  AE 
and  E  represent  the  energy  loss  in  the  AE  detector  and  in  the  E 
detector  respectively. 


AE  =  At 
E  t 


3. 4. 2.1 


where  At  and  t  are  the  thicknesses  of  the  respective  plastics. 


AE  =  a  AP 


3. 4. 2. 2 


E  =  b  P 


3. 4. 2. 3 


where  —  and  ~  are  the  gains  of  the  AE  detector  and  the  E  detector, 
a  b 

Substituting  these  equations  in  eq.  3. 4. 2.1  we  have: 


P_  At 
AP  t 


3. 4. 2.4 


The  absolute  gains  could  be  fixed  arbitrarily  as  — —  and  — —  for  AE 

6+1 

and  E,  or  equivalently  as  — and  (6+1)  respectively. 


ETOT 


6 

6+1 


AP  + 


1 

6+1 


P 


The  normalization  for  this  energy  was  determined  by  using  the 


ETOT  v/s  ®recoil  plot  for  (P*7*0)  and  lowing  the  energy  of  the 
bending  point  from  the  kinematics. 

The  empirical  relations  thus  arrived  at  are: 

ETOT  (LEFT)  -  0.45 -PE  +  0.22 -APE  3. 4. 2. 5 

ETOT  (RIGHT)  =  0.48 *PE  +  0.18 -APE  3. 4. 2.6 

3.4.3  Identification  of  (p,y)  Events 

The  off-line  analysis  was  done  using  a  DEC  VAX-11/780  computer 
with  a  VAX/VMS  operating  system.  The  analysis  program  used, 

’KIOWA*,  is  described  in  detail  in  Ste-81a.  The  program  was 
modified  and  adapted  for  the  VAX-11/780.  This  program  read  a  set  of 
control  cards  representing  the  software  cuts,  imposed  these  cuts  on 
the  data  read,  chose  the  appropriate  values  of  ADC,  TDC,  MWPC  co¬ 
ordinates  and  other  user  defined  parameters  based  upon  the  cuts,  and 
generated  one-dimensional  and  two-dimensional  plots.  In  addition, 
it  kept  track  of  user  defined  scalers  which  were  essential  in 
evaluating  the  cross  section.  User  defined  scalers  included  (p,  y) 
events  for  three  angles,  split  into  contributions  from  spin  up,  spin 
down  and  spin  off. 

Fig.  3.3  gives  the  flowchart  underlying  the  analysis  program. 

By  convention,  side  was  determined  by  a’s.  LI  refers  to  the 
LEFT  1  event,  which  in  turn,  refers  to  the  most  forward  y  on  the 
RIGHT  side.  L2  refers  to  the  middle-angle  y,  and  L3  to  the  most 
backward  y  on  the  RIGHT  side. 

The  pulser  events  were  first  scaled  and  then  eliminated  from 
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the  main  course  of  analysis.  The  rest  of  the  events  were  separated 
into  spin  up,  spin  down,  and  spin  off,  and  scaled.  The  next  check 
was  to  determine  whether  it  was  a  LEFT  or  RIGHT  event.  If  the  event 
belonged  to  neither,  it  was  discarded.  The  analysis  of  LEFT  and 
RIGHT  took  the  following  identical  steps. 

The  events  labelled  'MM'  (described  earlier)  were  scaled  and 
eliminated.  A  check  was  made  to  see  whether  the  event  was  LI,  L2  or 
L3  (Rl,  R2,  R3  for  RIGHT  side).  Once  this  was  established,  the 
following  identical  steps  were  taken  for  LI,  L2,  and  L3  (Rl,  R2, 

R3). 

The  scatterplot  of  AE  vs  ETOT  was  formed  which  helped  separate 
charge-1  particles  from  charge-2  particles.  The  clear  separation 
achieved  between  protons  and  a's  may  be  seen  in  Fig.  3.4.  The 
unwanted  protons  were  eliminated  with  AE-ETOT  cut.  Fig.  3.5  gives 
the  spectrum  after  eliminating  the  protons. 

The  next  step  was  to  generate  the  scatterplot  of  ETOT  vs  time- 
of-f light  for  the  a-particles  (Fig.  3.6).  Of  the  two  distinct 
peaks,  the  prominant  one  represented  the  REALS  and  the  other  the 
ACCIDENTALS.  The  term  ACCIDENTALS  refers  to  the  accidental 
coincidence  between  a(y)  from  a  beam  burst  and  the  y(  a)  from  another 
beam  burst.  The  ACCIDENTAL  contribution,  when  subtracted  from  the 
REALS  gave  the  TRUE  number  of  coincidences  in  one  beam  burst. 

The  next  step  was  to  generate  the  Cherenkov  pulse  height 
histogram  corresponding  to  REALS  (Fig.  3.7)  and  ACCIDENTALS  (Fig. 
3.8.  These  two  histograms  helped  elucidate  the  fact  that  the 
prominant  peak  in  the  beginning  channels  of  the  Cherenkov  pulse 
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height  spectrum  was  accidentals  due  to  general  background  from  the 
titanium  and  copper.  This  was  cross-checked  by  obtaining  the 
corresponding  spectra  for  the  titanium  run  (figs.  3.8a  and  3.8b). 
Consequently,  the  counts  from  the  initial  channels  of  the  Cherenkov 
pulse  height  histogram  were  eliminated  so  as  to  obtain  clear  spectra 
of  ETOT  vs  scattering  angle.  This  in  turn  necessitated  re¬ 
evaluating  the  Cherenkov  counter  efficiency  because  of  the  change  in 
the  threshold.  This  was  done  by  determining  the  change  in  the 
efficiency  due  to  the  software  cut  on  the  Cherenkov  PH  histogram  for 
d(p>Y)3He  which  we  did  along  with  the  tritium  runs  (Appendix  D). 

Fig.  3.9  gives  the  Cherenkov  PH  histogram  in  the  case  of  deuterium 
at  =  90°.  The  cross  section  in  this  case  is  larger  and  the 
counts  are  more  numerous,  which  helped  evaluate  the  efficiency 
factor  with  good  statistical  precision.  Note  also  that  this 
histogram  was  obtained  with  the  (p,y)  cut  in  the  energy  vs 
scattering  angle  histogram  and  therefore  the  events  were  due  to  y's 
only.  The  vertical  line  in  the  histogram  represents  the  pulse 
height  cut  made  for  the  tritium  case.  If  N1  represents  the  events 
to  the  left  of  the  line  and  N  the  total  events,  the  efficiency 
factor  involved,  n  is  given  by: 


3.4. 3. I 


The  final  step  in  the  analysis  was  to  generate  the  locus 
spectrum  (ETOT  vs  recoil  angle).  Figs.  3.10  and  3.11  represent  the 
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locus  spectrum.  Clear  separation  between  (p,y)  events  and  (p,7r°) 
events  may  be  noticed.  A  kinematic  check  was  made  to  reaffirm  the 
(p,y)  peak  obtained.  The  Cherenkov  counter  angular  acceptance  was 
calculated  from  geometry.  The  corresponding  angular  region  for  a’s 
was  found  from  the  kinematical  correspondence.  Again  from 
kinematics,  we  knew  the  energy  region  corresponding  to  this  peak. 
Since  the  energy  was  calibrated  in  absolute  terms  as  described 
earlier,  we  could  distinguish  the  (p,y)  peak  from  the  energy  point 
of  view  as  well.  The  (p,ir°)  locus  provided  additional  consistency 
checks.  The  kinematical  plot  (fig.  3.12)  illustrates  how  the  (p ,  tt°  ) 
locus  provided  the  consistency  check. 

3.5  THE  CROSS-SECTION  EVALUATION 

The  data  were  analyzed  for  the  TiT  runs  as  well  as  for  the 
blank  titanium  runs.  Tables  3.5a  and  3.5b  give  typical  summaries 
after  analyzing  a  number  of  runs  and  adding  up  the  appropriate 
counts  with  proper  normalization  factors.  The  program  'XSECT'  was 
written  to  evaluate  the  cross  section  and  the  analyzing  power.  The 
following  summarizes  the  way  in  which  the  cross  section  was 

evaluated. 

3.5.1  The  (p,Y)  Events 

The  total  number  of  events  corresponding  to  the  (p,y)  reaction 
were  calculated  as  described  earlier.  Let  that  be  denoted  by  'A*. 
Note  that  ’A'  was  the  net  contribution.  That  is,  the  accidental 
contribution  was  subtracted  from  the  reals  in  arriving  at  ’A*. 


Note 
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TABLE  3.5a 

OFF-LINE  ANALYSIS  OUTPUT  -  POLARIMETER  COUNTS 


POLARIMETER 

COUNTS 

TARGET 

LEFT 

RIGHT 

SPIN  UP 

SPIN  DOWN 

SPIN  OFF 

SPIN  UP 

SPIN  DOWN 

SPIN  OFF 

TiT 

73963248 

41873708 

24461808 

41350880 

73617136 

24307816 

(1928551) 

(950999) 

(633684) 

(934497) 

(1819933) 

(611107) 

Ti 

29577746 

17046344 

24041232 

16587889 

29870354 

23867232 

(681173) 

(342670) 

(767857) 

(328604) 

(646987) 

(735271) 

No.  of  incident  protons  (TiT)  =  2.27  x  1015  (up);  2.28  x  10 15  (down); 
0.96  x  10 15  (off). 

No.  of  incident  protons  (Ti)  =  0.91  x  10 15  (up);  0.93  x  10 15  (down); 
0.94  x  1015  (off). 

Note:  Numbers  in  parenthesis  refer  to  accidentals. 

CM 

Angles  covered  by  the  run  are  0  =  34°,  78°,  99°. 

Y 
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MWPC  efficiency  (Ti)  =  69%  (LEFT);  68%  (RIGHT). 


OFF-LINE  ANALYSIS  OUTPUT  -  (p,y)  EVENTS 
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also  that  'A'  refers  to  the  counts  in  the  TiT  run.  This  means  that 
background  from  Ti  was  present  in  'A'.  Later  on,  we  will  see  how 

that  contribution  was  subtracted  looking  at  the  data  for  a  Ti  run. 

3.5.2  The  Total  Beam  Current 

The  total  number  of  protons  incident  on  the  target,  which 
generated  ’A’  events,  were  taken  from  the  polarimeter  data  in  the 
trailer  block.  Eq.  2.8.1  was  used  to  calculate  the  incident  proton 

count  rate. 

3.5.3  The  Target  Thickness 

The  target  thicness  was  evaluated  as  described  in  Appendix  A. 
The  thickness  of  Ti  and  Cu  was  the  same  for  the  TiT  and  Ti  targets. 
The  tritium  thickness  is  denoted  by  T2  and  the  Ti-Cu  thickness  by 

Tr 

3.5.4  The  Solid  Angle 

The  solid  angle  was  determined  by  the  Cherenkov  counter,  since 
the  relative  distances  of  the  wire  planes  and  the  Cherenkov  counters 
were  adjusted  such  that  the  Cherenkov  counter  image  was  always 
inside  the  wire  planes.  Let  'C'  denote  the  solid  angle,  A^  the 
front  face  area  of  the  counter  and  d^  its  distance  from  the  target. 


3.5.5  The  Wire  Chamber  Efficiency 


Section  3.3  describes  the  significance  of  the  wire  chamber 
efficiency.  Since  there  were  two  planes  (x  and  y)  per  wire  chamber, 
an  efficiency  was  calculated  for  the  x  and  y  planes.  Let  the 
effective  efficiency  be  represented  by  'D'. 

U  -  events  without  'XMM*  «AND  •  events  without  'YMM' 

Total  events 

where  'XMM"  refers  to  the  x  plane  'multiple'  or  'miss'  and  'YMM' 
refers  to  the  y  plane  multiple  or  miss. 

3.5.6  The  Cherenkov  Counter  Efficiency 

The  Cherenkov  Counter  efficiency  involved  two  factors:  one  is 
the  one  described  in  Appendix  C,  which  was  a  function  of  y  energy 
and  the  distance  of  the  Cherenkov  counter  from  the  target.  Let 
this  factor  be  denoted  by  E.  The  other  one  was  a  function  of  the 
software  cut  applied  to  the  Cherenkov  counter  pulse  height 
histogram  (section  3.4.3).  Let  this  factor  be  denoted  by  F. 

3.5.7  The  System  Live-Time 

Section  2.10  describes  the  significance  of  system  live-time.  A 
major  portion  of  the  system  dead-time  was  contributed  by  the 
computer  processing  the  event.  The  rest  of  it  was  contributed  by 
the  LATCH  attached  to  the  EVENT  FAST.  The  overall  dead-time  was 
determined  making  use  of  a  pulser.  The  pulser  was  fed  to  drive  the 
LED's  attached  to  the  plastic  scintillators  and  the  Cherenkov 
counter.  The  pulser  pulses  were  also  fed  directly,  along  with  the 


'E'  signal,  to  an  'OR'  gate  whose  output  formed  part  of  the  input  to 
the  ’AND'  gate  which  determined  'LEFT'  or  'RIGHT'  'DCR'.  The  pulser 
feed  was  made  proportional  to  the  beam  current  by  setting 
coincidence  between  the  pulser  generator  output  and  the  polarimeter 
'L+R'  output,  Pulser  pulses  presented  to  the  system  were  scaled  and 
recorded  in  the  trailer  block.  The  pulser  event  simulated  a  good 
event,  thus  getting  recorded  as  if  it  was  a  good  event.  However 
pulser  'DCR'  demarcated  them  as  pulser  events.  A  typical  value  of  a 
pulser  present  was  0.16  M  while  a  total  number  of  1.12  M  events  were 
recorded,  of  which  3601  represented  (p,y)  events  of  0^’m-  =  54°  for 
both  sides.  The  overall  system  dead-time  was  a  function  of  the 
pulser  accepted  by  the  system  ( 'DCR'  test)  and  the  pulser  presented 
to  the  system  (scaled).  If  G  denotes  the  system  live-time,  we  have: 

G  =  pulser  accepted 
pulser  presented 

3.5.8  The  'TI'  Run  Parameters 

Let  Aa  denote  the  number  of  events  falling  in  the  same  software 
region  as  for  (p,y)  events  due  to  tritium.  The  number  of  incident 
protons  are  represented  by  B  the  target  thickness  by  T  A,  the  solid 
angle  (same  as  the  one  for  the  tritium  run)  by  C,  the  wire  chamber 
efficiency  by  Di,  the  Cherenkov  efficiency  factors  by  E  and  F  (same 
as  the  ones  for  the  tritium  run),  and  the  system  live-time  by  G^ 


3,5.9  Extraction  of  Events  Due  to  Tritium  and  Cross  Section 
Calculation 


Of  the  total  number  of  events  A,  not  all  stand  for  events  due 
to  tritium.  Letting  A'  denote  the  number  of  events  due  to  Ti-Cu,  we 
have: 


Af  =  Aj 


3. 5. 9.1 


The  number  of  events  due  to  tritium  alone  is  denoted 


by  at 


At  =  A  -  A* 


3. 5. 9.2 


The  differential  cross-section  x  is  given  by: 


x  =  — g'-.~  Y~'~  »  c  ♦  D  ~»~E  "  F~~T~G~  barns/steradian  3. 5. 9. 3 

where  T 2'  is  the  number  of  tritium  nuclei/cm2  of  the  target. 

3.6  THE  STATISTICAL  UNCERTAINTY  IN  THE  CROSS-SECTION  VALUE 

The  standard  deviation  in  the  value  of  x  may  be  obtained  from 
eq •  3. 5. 9. 3. 


a  =  x 
x 


at2  b 


1/2 


3.6.1 


where  D  =  —  and  G 
e 


=  & 
h  * 


3.7  THE  ANALYZING  POWER  EVALUATION 


When  we  followed  the  method  described  above  after  splitting  the 
total  events  into  contributions  due  to  spin  up  and  spin  down,  we  got 
the  spin-dependent  cross  sections.  This  was  done  for  both  the 
'LEFT*  side  and  ’RIGHT'  side.  The  analyzing  power  was  evaluated 
from  these  cross  sections.  Let  AP  represent  the  analyzing  power. 


AP 


L 

P 


h±Z\±lhil^M 


3.7.1 


where  X^  etc.  are  differential  cross-sections  for  LEFT,  SPIN  UP 
etc.,  and  P  is  the  average  beam  polarization  given  by: 


P 


VU  +  I^D 

XU  +  ID 


3.7.2 


where  P^(P^)  is  the  beam  up  (down)  polarization  given  by  eq.  2.8.2. 
Iy(ID)  are  the  beam  intensities  for  up  (down)  polarization.  The 
variance  in  the  analyzing  power  is  given  by: 

a f  2  =  ~~  (AP)  2  a_2  +  — - - - 

P2  P  P  2(m+n+r+s)  ** 

{(n+s)2  (a  Z+a  2)  +  (m+r)  2  (a  2)  }  3.7.3 

my  ns 

where  m  =  a^,  n=  a^,  r  =  0^,  s  =  °kf  (°Lf  etc.  refer  to  the 

standard  deviation  in  XT  .  etc.)  and 

L  t 
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a  2  =  - -  (I  2  a  2  +  I  2  a  2) 

P  Uu+ID>2  «  PU  D  PD 

( p  _p  )  2 

+ - S — (I  2  a  2  +  I  2  a  2)  3. 7. A 

(Iu+Iq)4  D  IU  U  ID 

substituting  eq.  3.7.4  in  eq.  3.7.3  we  obtained  the  standard 
deviation  in  the  analyzing  power. 
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CHAPTER  4 

RESULTS  AND  CONCLUSIONS 

4.  1  THE  RESULTS 

The  experimental  results  for  the  differential  cross  section  and 
the  analyzing  power  are  tabulated  in  table  4. 1  The  errors  quoted 
represent  only  the  statistical  uncertainty.  These  are  in  the  range 
from  6%  to  40%  for  cross  section  and  9%  to  48%  for  analyzing  power. 

The  results  were  converted  to  cross  section  values  for  the 
inverse  reaction,  using  the  principle  of  detailed  balance,  at  the 
corresponding  incident  y  energy  determined  by  the  conservation  of 
total  CM  energy.  Table  4.2  gives  the  results  for  the  inverse 
reaction.  This  conversion  was  done  so  that  the  results  could  be 
compared  with  the  existing  inverse  reaction  cross  section  values 
(Arg-75,  Pic-70,  Are-79,  Sch-83,  Kie-73).  These  correspond  to  data 
from  SACLAY,  FRASCATI,  BONN,  MIT,  and  ILLINOIS  respectively. 

Results  were  also  compared  with  forward  reaction  cross  section 
values  from  Did-70  (ORSAY). 

The  results  are  plotted  in  Fig.  4.1  to  4.3a.  Fig.  4.1  gives 
the  angular  distribution  of  the  cross  section  at  y-energy 
T_^  =  246  MeV.  The  results  are  shown  along  with  other  data  for 
comparison.  Figs.  4.2,  4.2a,  and  4.2b  give  the  excitation  function 
including  the  data  from  six  other  groups.  Fig.  4.3  gives  the 
analyzing  power  values  of  T^  =  300  MeV. 

There  is  fairly  good  agreement  between  the  results  of  the 
present  experiment  and  the  results  of  other  groups  with  the 
exception  of  those  of  Frascati  and  Illinois.  This  agreement 
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TABLE  4.1 

EXPERIMENTAL  RESULTS  FOR  CROSS  SECTION  AND  ANALYZING  POWER 


INCIDENT 

PROTON  ENERGY 
(MEV) 

eCM 

Y 

,  CM 

(nb/sr) 
a  M 

ANALYZING  POWER 

227 

90 

31.91  ±  2.21 

0.1167  ±  0.0888 

300 

54 

64.62  ±  3.78 

0.1050  ±  0.0971 

300 

65 

49.93  ±  3.10 

0.1180  ±  0.0888 

300 

78 

22.36  ±  1.69 

0.4355  ±  0.1301 

300 

90 

16.32  ±  1.77 

0.0337  ±  0.0305 

300 

99 

7.69  ±  1.50 

-0.6272  ±  0.3010 

300 

118 

4.94  ±  1.98 

“0.5842  ±0.3210 

375 

60 

37.69  ±  4.42 

0.4090  ±  0.2222 

375 

91 

14.34  ±  3.23 

-0.6147  ±0.3564 

126 


TABLE  4.2 

CROSS  SECTION  VALUES  CONVERTED  USING  DETAILED  BALANCE 


INCIDENT 

GAMMA  ENERGY 
(MEV) 

eCM 

p 

.  CM 

fa  (nb/sr) 

191 

90 

474.8  ±  32.8 

24b 

54 

786.8  ±  46.0 

246 

65 

608.0  ±  37.8 

24b 

78 

272.2  ±  20.5 

1  246 

90 

198.7  ±  21.6 

246 

99 

83.6  ±  18.3 

246 

118 

60.2  ±  24.1 

302 

60 

395.1  ±  46.2 

302 

91 

150.3  ±  33.8 
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indicates  that  time  reversal  is  not  violated  in  the  electromagnetic 
interaction  of  hadrons.  This  aspect  is  discussed  later. 

Figs.  4.1,  4.2a  and  4.2b  also  give  theoretical  curves  for  the 
cross  section  values,  calculated  in  the  Distorted  Wave  Impulse 
Approximation  (DWIA).  This  is  discussed  in  the  next  section. 

4.2  DISTORTED  WAVE  IMPULSE  APPROXIMATION  (DWIA) 

The  calculations  performed  are  based  upon  the  quasi-deuteron 
model  first  proposed  by  Levinger  (Lev-51).  This  model  considered 
the  total  nuclear  wavefunction  as  the  product  of  the  wavefunction 
for  two  nucleons  very  close  together  and  the  wavefunction  for  the 
remaining  nucleons.  The  remaining  nucleons  provided  a  potential  in 
which  the  deuteron  moves  and  they  scatter  the  high  energy  nucleons 
of  the  deuteron.  The  wavefunction  for  the  A-nucleon  system  was 
expressed  as  a  product  of  three  terms:  the  first  term  representing 
the  motion  of  center  of  mass  of  the  quasi-deuteron,  the  second  term 
representing  the  wavefunction  of  the  quasideuteron,  and  the  third 
terra  representing  the  wavef unctions  of  the  A- 2  nucleus. 

¥  (l,  2,  ...A)  =  exp  (i£f  •  r'  )  1)^  (r  )  $(3,  ...A)  (4.1) 

where  'F  is  the  total  wavefunction,  exp(ik’  *r ’ )  stands  for  the  motion 
of  the  center  of  mass  of  the  deuteron,  ^(r)  stands  for  the 
quasideuteron  wavefunction  and  c|>  stands  for  the  wavefunction  of  the 
remaining  A-2  nucleons.  The  wavenumber  for  the  relative  motion  of 
proton  and  neutron  is  represented  by  k  and  the  internucleon  distance 
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is  represented  by  r.  The  cross  section  for  the  photonuclear  process 
was  obtained  in  terms  of  the  cross-section  for  the 
photodisintegration  of  the  quasideuteron,  given  the  velocity  v^  of 
the  quasideuteron  in  the  direction  of  observation  and  v^  in  the 
perpendicular  direction.  A  proton  observed  at  an  energy  T  outside 
the  nucleus  is  scattered  at  an  energy  T*  by  the  quasideuteron.  The 
same  proton,  observed  at  an  angle  0  is  scattered  at  an  angle  0'  by 
the  quasideuteron.  The  differential  cross  section  as  a  function  of 
T,  was  given  by: 


do  (T)  =  /do(T'  )  P(vx)  dvx  (4.2) 

where  P(v^)  is  the  distribution  function  for  the  velocity  component 

v  . 

x 

The  differential  cross  section  as  a  function  of  angle  was  given 
by: 


d o(  6)  =  /do(  0’)  P(vy)  dVy  (4,3) 

where  P(v^)  is  the  distribution  function  for  the  velocity  component 

v  • 

y 

This  method,  with  modifications,  was  applied  to  calculate  the 
cross  section  for  the  reaction  pd  -►  frr+  by  Ruderman  (Rud-52).  In 
this  approach  the  cross  section  was  calculated  in  terms  of  the 
cross  section  for  the  primary  reaction  PP  d  The  neutron  in  the 
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deuterium  acts  as  a  spectator.  The  triton  wavefunction  was 
expressed  in  terras  of  the  deuteron  wavefunction  on  the  assumption 
that  the  neutron  and  proton  that  participated  directly  in  the 
production  of  the  triton  are  very  close.  Thus  the  following 
expression  for  the  cross  section  (Blu-54)  was  obtained: 


t  ir 


=  [f< 9)  ]2  * e~?V [Vpp  *  d  U (pp  +  d,,+ j] 

Pd  tt  t  d  -1 


(4.4) 


where  v’s  are  the  relative  velocities  and  E’s  the  total  energies  in 
the  CM  frame.  The  quantities  inside  the  square  brackets  were 
evaluated  in  the  CM  system  of  the  two  protons.  The  energy  of  the 
proton  was  determined  for  producing  a  meson  of  momentum  cf.  The 
quantity  f ( 0)  is  the  form  factor,  given  by: 


f  (  9) 


/%(*)  exp  (i  -^2-) 
n 


^(0,x) 

Vo> 


dx 


(4.5) 


where  ^  and  are  the  deuteron  and  triton  wavef unctions ,  x  the 

■*  1  -►  1  -»■ 

neutron-proton  relative  co-ordinate  and  A  =  -j  k  -  -j  q.  The  momentum 
of  the  incident  proton  is  given  by  It. 

This  model  was  modified  by  Fearing  to  incorporate  distortion  in 
the  initial  and  final  channels  (Fea-74,  Fea-75a)  and  to  include  a 
better  treatment  of  wavef unctions.  The  model  was  extended  to 
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include  the  reaction  Pd  *►  %e y  (Fea-75b),  with  the  difference  that 
there  is  no  outgoing  wave  distortion.  Described  below  are  the  main 
features  of  the  model  and  the  calculations  done  in  the  present  work 
to  get  the  cross  section  for  P  +  T  +  y  +  using  DWIA. 

4.3  SALIENT  FEATURES  OF  FEARING 'S  DWIA  MODEL 

The  cross  section  for  the  reaction  P  +  T  ■*  ^e  +  y  could  be 
expressed  as  a  product  of  the  cross  section  for  the  sub-process 
P+n  d  +  y,  a  kinematic  factor,  a  spin-isospin-antisymmetrisation 
factor,  and  a  form  factor. 

P  +  T  ■*  y  +  a  P  +  n*y+d 

)  -  KF.g.|F(A)|2  (||)  (4.6) 

PTcm  Pn  cm 

where  KF  is  the  kinematical  factor,  g  the  spin  factor,  and  F(A)  the 
form  factor. 


rd  a 


KF 


E 

q 


E 

a 


(4.7) 


where  E's  and  P’s  are  the  energy  and  momentum  in  the  PT  cm  frame. 
Eq .  refers  to  the  target  neutron: 


Eq  =  (q  2  +  m^  2 


i/2 

) 


(4.8) 


E  ,  =  E  +  E  -  E 

d  P  q  y 


(4.9) 
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fq'  is  the  effective  momentum  involved,  determined  by  taking 
into  account  the  fact  that  the  sub-process  takes  place  off-mass- 
shell. 

The  spin  factor  is  given  by: 


g 


Z 

spins 


Z 

|  3  3  3  <3 

I  d  o  i  d 


T  | 
01 


3  3  ><  a 

0  1  a 


3><33  | 

d  0  1 


a  a  > 
P  T 


2 


z  |<  ed 

spins 


I  e  e  >  I2 

0  1 


(4.10) 


where  3’s  refer  to  the  spin-isospin  wavefunct ions  of  those  involved 
in  the  process  pn  dy,  0  refers  to  the  incoming  proton,  1  to  the 
neutron,  and  a’s  refer  to  spin-isospin  wavefunct ions  of  those 
involved  in  the  process  PT  ya.  Both  the  numerator  and  the 
denominator  are  calculated  by  summing  over  the  final  spins  and 
averaging  over  the  initial  spins.  TQ1  is  the  two-body  T-matrix. 

The  value  of  g  for  the  reaction  P  +  T  +  y  +  ^lle,  was  calculated  and 
found  to  be  equal  to  1/18. 

The  form  factor  is  given  by: 


F(  A) 


2  ,  3  *  3  »  3  3  3  iA*v 

—  /  dydydydy  6  (y  +y)e 

/3  0  1  2  3  2  3 


(4.11) 
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where  unpriraed  and  primed  quantities  stand  for  those  before  and  after 

the  collision.  If  rQ,  r  p  r  2,  r3  represent  the  co-ordinates  of  the 

initial  four  nucleons  and  r^,  r r^,  r  ^  those  of  the  final  four 

nucleons,  then  y.  =  r.  -  R  and  y!  =  r!  -  R  where  R  is  the  center 
i  1  s  Ji  i  s  s 

1  PT 

of  mass  of  the  spectator  nucleons  n  and  p.  A  =  -=■  (P  -  - P  ) 

Zm  l  j  *y  Jl  i.  cm 

is  known  as  the  effective  momentum  transfer,  x  q  =  r  q  -  R^  where  RT 
is  the  triton  c.m.  co-ordinate.  Vp^  is  the  relative  velocity  between 
proton  and  triton  and  v^  that  between  proton  and  neutron.  $'s  are 
antisymmetric  wavefunctions .  The  F’s  are  the  distortion  functions 
which  are  measures  of  the  deviation  of  the  scattering  waves  from 
plane  waves.  F  is  found  by  integrating  the  scattering  potential 
along  the  projectile  path.  The  potential,  in  turn,  is  expressed  in 
terms  of  the  nuclear  density  function,  the  forward  scattering 
amplitude  and  the  total  cross  section  for  p-p  elastic  scattering  and 
p-n  elastic  scattering. 


F(x,v) 


exp 


+  i)  +  2  o  ( n  +  i) 
n  n 


a 

/  p(x— vi)dT 
J  0 


(4.12) 


where  is  the  total  cross  section  for  p-p  scattering,  and  on  is  the 
total  cross  section  for  p-n  scattering.  n  is  the  ratio  of  the  real 
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to  imaginary  parts  of  the  forward  amplitude  for  p-p  scattering  and  n 
the  analogous  quantity  for  p-n  scattering.  The  charge  density 
function  is  represented  by  p.  It  is  given  by  the  following 
equation: 


p(r) 


(r-R)/a 

1  +  e 


(4.13) 


where  pQ  is  the  central  density,  ’a'  the  skin  thickness  and  ?R'  the 
radius  obtained  from  the  electron  scattering  data.  The  central 
density  is  given  by: 


P 

o 


_ 3 _ 

4 ir  (R3  +  tt^^O 


(4.14) 


4.4  THE  WAVEFUNCT IONS  OF  d,  t  AND  a 

The  form  factor  evaluation  requires  knowledge  of  the 
wavef unctions  of  the  deuteron,  triton  and  alpha  (eq.  4.11).  The 
Hulthen  wavefunction  (Fea-75a)  was  used  for  the  deuteron.  The 
Hulthen  deuteron  wavefunction  is  given  below. 


(4.15) 
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1/2 

where  N  =9.8  MeV 
d 

8  =  45  MeV 
Y  =  270  MeV. 

A  correlated  Gaussian  wave  function  (Fea-77)  was  used  for  the 
triton.  The  following  is  the  functional  form. 

4>t  (x0,xpx2)  =  Nt  e  a2u[l  -  C  exp(-82|x  0-x  J  2)  ] 

x  [l  -  C  exp(- $2|x  Q-x  2|  2)  ][l  ~  C  exp(-82|x  x-x  2|  2)  ]  (4.16) 

where  u  =  |xQ-x  J  2  +  |x  A-  x2|  2  +  |xQ-  x2|  2  , 

Nt  =  1.83  x  10 6  MeV3  , 
a  =  62.8  MeV  , 

8  =  232  MeV  ,  and 
C  =  0.925. 

A  correlated  exponential  wavefunction  was  also  used  for  triton.  The 
form  is  as  given  below:. 

4>3(x0,x  px2)  =  Nt  e  “*[1  -  C  exp(- 82|x  Q-x  J  2)  ] 

x  [l  -  C  exp(-  82|x  q-x  2|  2)  ][l  -  C  exp(-  82|x  ^x  2|  2)  ] 

I  -►  >|  |  ■>  "►!  (->■  -►! 

where  w  =  Ixq-Xj^I  +  |xA~x2|  +  |Xq-x2|. 

A  correlated  Gaussian  function  was  used  for  ^e,  which  is  given 


below: 
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(xQ,x  a,x2,x3)  =  N^e  32u[l  -  C  exp(-y2|x0-x  J  2)  ] 
x  [l  -  C  exp(- y2|x  q-x 2|  2)][l  -  C  exp(- y2|x  Q-x  3|  2)  ] 
x  [l  -  C  exp(-y2|x  A-x2|  2)  ][l  -  C  exp(-y2|x  A~x3|  2)  ]  (4.17) 

x  [l  -  C  exp(  y2|x2-x  3|  2)  ] 

where  =  9.55  x  101*  MeV  3,  3  =  166  MeV,  and  y  =  615  MeV. 

Taking  a  close  look  at  eqn.  4.11  it  may  be  seen  that  the  form 
factor  is  the  fourier  transform  with  respect  to  the  momentum  transfer 
of  the  overlap  of  the  triton  and  a  wavefunctions  with  the  deuteron 
wavef unction,  modified  by  the  PT  and  Pn  distortion  function.  This 
means  that  in  order  to  have  an  appreciable  cross  section,  there  has 
to  be  a  good  amount  of  overlap  between  the  intermediate  deuteron  and 
final  state  nuclear  wavefunctions,  a  condition  which  is  satisfied  for 
low-mass  nuclei. 

4.5  DATA  FOR  THE  REACTION  P  +  n  •»  y  +  d 

This  section  elaborates  on  the  last  terra  in  eq.  4.6,  namely  the 
experimental  cross  section  values  used  for  the  reaction 
P  +  n  -►  y  +  d.  Reaction  amplitudes  for  the  reaction  P  +  n  -►  y  +  d 
can  be  substituted  for  the  interaction  vertex,  thereby  eliminating 
the  ambiguities  associated  with  the  Single  Nucleon  Model  (SNM)  and 
the  Two  Nucleon  Model  (TNM).  In  the  SNM,  the  proton  emits  a  y  and 
then  is  captured  by  the  triton  to  form  ^e.  The  entire  momentum 


transfer  is  taken  up  by  a  single  nucleon,  thereby  involving  high 
momentum  components  of  the  wavef unction,  which  are  not  yet  understood 
properly.  In  the  TNM,  the  incoming  proton  shares  its  momentum  with 
one  nucleon  of  the  triton  either  in  the  process  of  quasideuteron 
formation  or  A  isobar  production.  Here  momentum  transfer  is  shared 
by  two  nucleons,  which  eliminates  the  need  for  high  momentum 
components. 

The  experimental  cross  section  value  for  y  +  d  +  P  +  n  has  been 
measured  recently  with  high  precision  over  a  wide  energy  range  by 
four  groups  (Cra-82,  Are-84,  Bab-83,  Dou-76).  These  data  come  from 
MIT,  BONN,  TOKYO  and  LUND  respectively.  Taken  as  a  whole,  these  data 
sets  cover  an  energy  range  of  =  100  MeV  to  832  MeV  over  an  angular 

nvi 

range  of  0^  =15°  to  160°.  The  values  from  LUND  below  300  MeV  have 

been  discarded. 

Fig.  4.4  gives  a  plot  of  these  data  as  angular  distributions  of 
the  cross  section  for  various  incident  y  energies.  Plotted  also  are 
the  fitted  curves  at  appropriate  energies.  Fitting  was  done  taking  s 
cros  section  values  from  all  the  groups.  The  cross  section  values 
at  a  particular  energy  were  fitted  to  a  polynomial  of  the  form  given 
below: 


f§  =  a0  +  ai  cos 


PCM 


a  2  cos 


20 


PCM 


+  a  o  cos  J0 


PCM 


(4.18) 


where  a's  are  the  parameters  to  be  determined  and  0p^  is  the  proton 


scattering  angle  in  the  CM  frame.  Fitting  the  data  to  the  polynomial 
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150 


6  (CM) 
P 
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(4.18)  was  done  for  a  set  of  energies  and  the  values  of  'a'  were 
generated  as  a  function  of  E^.  Table  4.3  gives  the  coefficients  of 
the  polynomial  thus  generated. 

These  coefficients  were  then  fitted  as  a  function  of  energy. 
The  following  functional  forms  were  used. 


a  n  *  cOl  exp(  -c02*x)  + 


c03  -  c04  «x 

1  +  fx-c05,i5 
l-c06"J 


(4.19) 


ell  -  c  1 2  *x  +  c  1 3  *x  2 


(4.20) 


2  oi  oo  c23 

—  =  c2 1  -  c22  -x — 


(4.21) 


3  oi  oo  c33 
—  =  c31  -  c32  -x - 


(4.22) 


Table  4.4  gives  the  values  of  the  above  parameters.  The  C 
values  were  substituted  in  eqs.  4.19  to  4.22  to  generate  a's  which 
were  then  used  in  eq.  4.18  to  give  cross  section  at  any  energy  and 
angle.  This  cross  section  value  was  then  used  in  eq.  4.6  to  yield  a 
cross  section  value  for  P  +  T  -*•  ^He  +  y. 

4.6  lM| 2  VERSUS  t  REPRESENTATION  OF  THE  RESULT 


A  better  way  of  looking  at  the  result  than  the  versus  0 


300 

325 

350 

375 

400 

425 

450 

475 

500 

525 

550 

575 

600 

625 


152 


TABLE  4.3 

COEFFICIENTS  OF  THE  POLYNOMIAL 


ao 

al 

a  2 

a  3 

5.30  ±  0.07 

0.52  ±  0.32 

-1.61  ±  0.27 

0.90  ±  0.60 

4.34  ±  0.04 

0.83  ±  0.29 

-1.43  ±  0.25 

-02.0  ±0.60 

3.24  ±  0.04 

0.75  ±  0.20 

-0.94  ±  0.17 

-0.20  ±  0.40 

2.40  ±  0.05 

0.61  ±  0.14 

-0.94  ±  0.12 

-0.70  ±  0.30 

1.77  ±  0.03 

0.35  ±  0.12 

-0.38  ±  0.12 

-0.30  ±  0.20 

1.45  ±  0.02 

0.36  ±  0.11 

-0.38  ±  0.10 

-0.20  ±  0.40 

1.16  ±  0.02 

0.30  ±  0.11 

-0.26  ±  0.10 

-0.10  ±  0.20 

1.08  ±  0.02 

0.18  ±  0.18 

-0.41  ±  0.15 

0.30  ±  0.40 

0.93  ±  0.04 

0.13  ±  0.27 

-0.20  ±  0.34 

0.40  ±  0.70 

0.80  ±  0.03 

0.33  ±  0.26 

-0.09  ±0.33 

-0.10  ±  0.70 

0.72  ±  0.02 

0.14  ±  0.17 

0.07  ±  0.49 

0.40  ±  0.70 

0.69  ±  0.02 

0.11  ±  0.13 

-0.43  ±  0.39 

-0.10  ±  0.60 

0.60  ±  0.02 

0.34  ±  0.12 

0.19  ±  0.62 

0.30  ±  0.70 

0.51  ±  0.02 

0.60  ±  0.15 

1.08  ±  0.74 

0.90  ±  0.80 
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TABLE  4.4 

'C'  PARAMETERS  OF  EQUATIONS  4.17-4.20 


CU1 

C02 

C03 

C04 

C05 

C06 

132. 25± 

10.50 

0.39  ± 

0.01 

6.23  ± 

0.07 

0.0040  ± 

0.0002 

247.25  ± 

1.18 

14853  ± 

480 

Cll 

C12 

Cl  3 

0.6650  ±  0.0251 

0.0029  ±  0.0002 

4. 14x10” 6  ±  4. 10x10“ 7 

;  C21 

C22 

C23 

0.150  ±  0.001 

-0.1210  ±  0.0002 

-0.1210  ±  0.0002  j 

C31 

C32 

C33 

-0.0943  ±  0.0038 

-0.0002  ±  0.001 

5.9048  ±  2.3613 
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representation  is  |m|  2  versus  t  representation,  since  the  phase  space 
factor  is  removed  in  the  latter  representation.  |m|  2  is  known  as  the 
invariant  matrix  element  squared  (Wil-71)  which,  for  the  reaction 
Y  p  +  t  is: 

l«|2=  ^  ^2Sy  +  1)(2S„+  1)  (4.23) 

where  S^,  are  the  photon  and  a  spin,  S  the  square  of  C.M.  energy, 

P  ,  Pp,  are  in  the  C.M.  frame  of  reference.  't'  is  the  square  o^ 
Y  ^  cl  u  e 

the  4-moraentum  transfer.  Fig.  4.5  gives  such  a  representation  of  tha 
result.  It  may  be  seen  that  60°CM  and  90°CM  results  do  not  fall  on 
universal  curve  unlike  the  case  of  y  +  %e  -*■  P  +  d  (Cou-83). 

4.7  CONCLUSIONS 

The  angular  distribution  of  the  cross  section  at  246  MeV 
incident  energy  is  plotted  in  Fig.  4.1  along  with  the  data  from 
BONN,  MIT  and  SACLAY  and  the  present  work.  The  theoretical  curves 
calculated  from  DWIA  are  also  plotted.  Curve  A  and  curve  B  stand 
for  the  Correlated  Gaussian  wavefunction  used  for  triton,  A  standing 
for  the  case  without  distortion  and  B  with  distortion.  As  may  be 
seen,  the  distortion  decreases  the  cross  section  by  around  29%. 

This  means  that  the  distortion  effects  are  of  considerable 
importance.  Since  the  distortion  represents  dominantly  the 
absorption  of  the  incoming  proton,  the  cross  section  is  reduced  as 
expected,  when  distortion  is  included. 
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All  the  curves  reproduce  the  shape  very  well  except  at  backward 
angles,  which  is  characteristic  of  the  model  in  situations  like 
p+d  ■*  t+ir+  (Fea-75)  and  p+d  +  %e  +  y  (Cam-83).  This  could  be  because 
mechanisms  involving  more  than  two  nucleons  become  important 
resulting  in  the  reduced  momentum  transfer  which  results  in  the 
increased  form  factor  and  hence  the  increased  cross  section.  Since 
the  model  is  concerned  at  the  primary  mechanism  pn  dy,  thereby 
confining  itself  to  the  two-nucleon  mechanism,  it  is  natural  that 
the  theoretical  prediction  is  lower  than  the  experimental  value.  A 
typical  situation  is  at  140°  where  the  experimental  value  is  2.5 
times  that  of  the  theoretical  calculation. 

It  may  be  observed  that  the  distortion  introduces  a 
normalization  effect  without  changing  the  shape  of  the  curve.  Again 
this  is  characteristic  of  the  model,  observed  in  the  situation  of 
p+d  +  t+TT+  also.  This  is  understandable  enough,  since  the 
distortion  functions  are  evaluated  at  the  origin  and  so  distortion 
effect  should  be  independent  of  the  coordinates. 

Curves  C  and  D  represent  calculations  using  the  correlated 
exponential  wavefunction  for  the  triton,  C  representing  the  case  of 
no  distortion  and  D  the  case  with  distortion.  Curve  C  gives  very 
good  agreement  with  the  experimental  situation.  One  tends  to 
conclude  that  the  correlated  exponential  wavefunction  gives  the  best 
representation  of  the  triton. 

Fig.  4.2  gives  the  experimental  results  for  the  excitation 
function  at  60°  and  90°  C.M.  angles.  Plotted  also  are  the  data  from 
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BONN,  SACLAY,  MIT,  ILLINOIS  and  ORSAY.  Fig.  4.2a  gives  the  60° 
excitation  function  along  with  the  DWIA  calculation  results. 
Calculation  was  done  using  the  correlated  exponential  wavefunction 
for  the  triton.  It  may  be  seen  that  the  shape  of  the  data  is  well 
imitated  by  theory  except  at  large  energies  where  the  theory 
predicts  values  lower  by  a  factor  of  2  at  a  typical  value  of  E  ^  = 

300  MeV.  Again,  this  is  characteristic  of  the  theory  illustrated  in 
the  case  of  p+d  t+ir+.  The  multi-nucleon  mechanism  becomes 
important  at  high  energies  resulting  in  the  increased  form  factor 
and  the  cross  section.  Since  the  theory  does  not  take  this  into 
account,  it  is  falling  below  the  experimental  value  at  high 
energies. 

Fig.  4.2b  represents  the  excitation  function  at  90°  C.M.  The 
fall-off  at  higher  energies  is  sharper  than  compared  to  the  60° 
case.  Theory  underpredicts  the  cross  section  by  a  factor  of  3  at 
300  MeV,  compared  to  a  factor  of  of  2  at  60°  C.M. 

Fig.  4.3a  gives  the  experimental  value  of  the  excitation 
function  at  60°  C.M.  along  with  the  curve  representing  the 
calculation  in  Gari  and  He bach  theory.  Curve  A  represents 
calculation  including  the  A  effect.  Curve  B  represents  calculation 
excluding  the  A.  Both  the  curves  follow  the  shape  well  although 
predicted  values  are  lower  by  a  factor  of  1.7  at  200  MeV.  The  A 
effect  becomes  important  beyond  150  MeV.  The  theory  underpredicts 
the  cross-section  by  a  factor  of  1.3  at  200  MeV  when  the  A  effect  is 
included.  The  agreement  with  the  experiment  becomes  very  good 
beyond  225  MeV  when  the  A  is  included. 
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A  preliminary  discussion  of  Time  Reversal  Invariance  (TRI)  was 
given  in  Chapter  1.  We  have,  at  our  disposal  a  means  to  check  the 
TRI  by  checking  the  principle  of  Detailed  Balance  (eqn.  1.62).  As 
mentioned  in  Chapter  1,  the  comparison  of  the  magnitude  of  the  cross 
section  and  the  shape  of  the  angular  distribution  are  two  criteria  to 
check  TRI.  This  is  because  the  ratio  of  the  forward  cross  section  to 
the  backward  cross  section  is  given  by  a  quantity  independent  of  the 
scattering  angle.  Since  both  the  criteria  are  satisfied,  as  may  be 
seen  when  comparing  the  present  data  with  the  inverse  reaction  data, 
we  see  no  evidence  for  the  time  reversal  violation.  As  seen 
earlier,  y-N-N  vertex,  current  conservation  and  hermiticity  prevent 
and  T  violation  in  these  interactions  except  in  the  A  resonance 
region.  The  present  measurements  do  cover  the  A-resonance  region 
and  thus  do  amount  to  an  important  check  on  TRI.  A  similar 
conclusion  has  been  obtained  by  comparing  the  results  of  P  +  d  ■*  ^le 
+  y  to  the  inverse  reaction  values  (Cam-83).  The  results  for  the 
reaction  P  +  d  -►  %e  +  y  mentioned  in  Appendix  D  are  also  consistent 
with  the  above  claim. 

Fig.  4.3  presents  the  analyzing  power  of  the  reaction  p  +  T 

a  +  y  at  T  =  300  MeV  as  a  function  of  the  scattering  angle.  The 
P 

positive  values  at  small  angles  and  negative  values  at  large  angles 
are  reminiscent  of  the  analyzing  power  for  pd  +  %e y  (Cam-83).  The 
negative  value  of  the  analyzing  power  might  thus  be  attributed  to 
quasideuteron  contributions. 

Since  the  analyzing  power  involves  the  ratio  of  various 
amplitudes,  it  is  not  sensitive  to  the  overall  normalization,  unlike 


the  case  of  differential  cross-section. 

Fig.  4.6  presents  the  analyzing  power  for  the  reactions  p  +  T  ■+■ 
a  +  y  at  =  300  MeV,  compared  to  kinematically  shifted  analyzing 
powers  for  the  reactions  p  +  n  -►  d  +  y  and  p  +  d  %e  +  y.  The 
kinematic  correspondence  is  found  by  involking  the  Impulse 
Approximation  where  the  subprocess  amplitude  pn  dy  is  equated  to 
the  free  scattering  amplitude  pn  *  dy  by  equating  the  momentum 
transfers  as  well  as  the  centre  of  mass  total  energies  (Fea-75). 

The  magnitude  is  found  to  be  closest  at  forward  angles  to  p  +  d 
-*•  %e  +  y  and  at  backward  angles  to  p  +  n  d  +  y.  The  backward 
angles  results  however  also  agree  with  the  p  +  d  >  %e y  result  to 
within  one  standard  deviation.  It  thus  appears  that  the  Impulse 
Approximation  can  account  for  not  only  the  cross-section  but  also 


the  analyzing  power. 
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APPENDIX  A 


TARGET  THICKNESS  EVALUATION  FROM  P-T  ELASTIC  SCATTERING 

A  preliminary  description  of  the  target  construction  has 

already  been  given  in  Chapter  2. 

The  thickness  was  determined  by  elastically  scattering  protons 

from  the  TiT  and  Ti  targets  using  the  MRS  spectrometer  facility. 

The  known  cross-section  values  for  P+T+P+TatE  =  318  MeV 

P 

(Bli-82)  was  used  in  determining  the  tritium  content  of  the  TiT 
target.  Additional  measurements,  made  at  E^  =  500  MeV,  were 
compared  with  the  cross-section  values  extrapolated  from  the  results 
of  G.  Bizard  et^  al.  (Biz-80). 

A. 1  THE  MEDIUM  RESOLUTION  SPECTROMETER 

This  magnetic  spectrometer  consists  of  a  dipole  and  quadrupole. 
The  dipole  field  strength  is  such  that  the  maximum  momentum  of  the 
detected  particle  is  1.5  GeV/c.  A  dipole  with  a  60°  bend  angle 
follows  a  quadrupole  which  is  preceded  by  a  wire  chamber  and  a 
scattering  chamber,  looking  downstream.  The  wire  chamber  XqYq  is  of 
dimensions  5"  x  5”  with  64  wires  in  each  plane.  Beyond  the  exit  of 
the  dipole  come  two  wire  chambers,  each  of  dimensions  10"  x  40". 

The  X  plane  has  512  wires  and  the  Y  plane  has  128  wires.  The  bottom 
chamber  is  known  as  X^-Y 

comes  a  1"  thick  plastic  scintillator.  Wire  spacing  of  any  wire 
chamber  is  2  mm. 

The  large  angle  configuration  (LAC)  of  the  MRS  provides  for 


and  the  top  chamber  as  X^-Y^.  After  these 
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measurement  at  scattering  angles  of  20°  and  beyond.  A  small  angle 
configuration  (SAC)  is  achieved  by  pushing  the  quad  2  m  from  the 
target.  In  this  configuration,  measurements  can  be  made  from  20°  to 
the  left  to  12°  to  the  right.  The  solid  angle  acceptance  of  the  MRS 
for  angles  above  15°  is  2.5  msr.  The  corresponding  number  for 
angles  less  than  15°  is  1.5  msr.  The  combined  energy  resolution  of 
the  cyclotron,  the  beam  transport  system,  and  the  spectrometer  was 
measured  and  found  to  be  900  keV  FWHM  at  400  MeV  and  600  keV  FWHM  at 
200  MeV.  Details  of  the  MRS  may  be  found  in  (Ste-81,  Mil-83). 

A. 2  THE  EXPERIMENTAL  PROCEDURE 

The  experiment  was  done  at  the  T2  target  station  using  the  MRS. 
The  proton  beam  of  318  MeV  was  used.  The  protons  elastically 
scattered  from  titanium,  copper  and  tritium  were  analyzed  by  the 
MRS.  In  addition  we  had  to  contend  with  the  protons  elastically 
scattered  from  carbon  and  hydrogen  which  constituted  the  kapton 
window  for  the  target.  Note  that  the  target  enclosure  had  two 
layers  of  1  mil  kapton  window  on  either  side.  We  could  easily 
separate  protons  due  to  tritium  from  the  rest  of  the  protons  by 
choosing  the  appropriate  kinematical  region  (see  Fig.  A.l).  The 
MRS  was  set  at  LAC  mode  and  measurements  were  made  at  three 
scattering  angles,  03  =  18.04°,  22.01°  and  25.99°.  We  had  one  TiT 
run  and  one  Ti  run  for  03  =  18.04°,  five  TiT  runs  and  two  Ti  runs 
for  03^  =  22.01°,  three  TiT  runs  and  one  Ti  run  for  03^  =  25.99°. 

See  Table  A.l  for  a  run  summary. 
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A. 3  THE  MOMENTUM  ANALYSIS 

The  focal  plane  position  of  the  scattered  particle  is  a 

function  of  its  momentum.  This  meant  that  the  different  momentum 

groups  of  protons  could  be  identified  by  their  position  in  the  focal 

plane.  Focal  plane  position  is  given  in  terms  of  the  wire  numbers 

X  X  and  X^  hit  by  the  particle.  If  X  stands  for  the  focal 
O,  B  i  r 

plane  position,  we  have: 


359  •  XT  -  139  •  Xfi 
XF  =  ^  -  Xg  +  220 


0.07  (X  -  28) 2 
o 


(A.  1 ) 


A  kinematical  correction  depending  upon  the  value  of  Yq  needed 
to  be  applied  to  this  so  that  vertical  bands  corresponding  to 
different  momenta  could  be  obtained  while  plotting  Yq  versus  the 
modified  focal  plane  position  (denoted  by  X  ). 

r  K. 

XFK  =  Xf  +  BYo  (A. 2) 

where  B  depends  upon  the  scattered  proton  momentum.  Fig.  A. 2  gives 

a  plot  of  Y  versus  X„„  where  kinematical  correction  has  been  done 
o  IK 

for  triton.  Figures  A. 3,  A. 4  and  A. 5  give  the  one -dimensional  plots 

for  X  for  scattering  angles  18.04°,  22.01°  and  25.99°.  Fig.  A. 6 
r  K. 

is  the  corresponding  plot  for  the  Ti  run. 

A. 4  TARGET  THICKNESS  EVALUATION 
A . 4 . 1  The  Solid  Angle 


The  solid  angle  is  determined  by  the  window  put  on  the  X<j>-Y<j> 


P-T  ELASTIC  SCATTERING  DATA  SUMMARY 
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CHNLL  MEANS  LOW  CHANNEL  CORRESPONDING  TO  TRITIUM  PEAK. 
CHNLH  MEANS  HIGH  CHANNEL  CORRESPONDING  TO  TRITIUM  PEAK 
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Tj-  U0  00  GO  CVJTTUJCD  CD 

co  cn  cd  a>  r^uoin^J- 


450  498  546  594  642  690  738  786  834 
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S 


l  N  n  0  D 


400  464  528  592  656  720  784  848  912 

XFK  (mm) 
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plane.  The  window  was  such  that  50  mm  to  81  mm  of  the  X<J>  plane  and 
50  mm  to  86  mm  of  the  Y<f>  plane  contributed  to  the  solid  angle.  The 
X<}>  plane  was  situated  at  a  distance  of  85.09  cm  and  the  Y <|>  plane  at 
a  distance  of  87.63  cm  from  the  target.  From  geometry,  the  angle 
subtended  by  the  relevant  portions  of  the  Y<J>  plane  is  found  to  be 
15.4045  mr.  The  corresponding  quantity  for  the  XcJ>  plane  is  25.6705 
mr.  These  numbers  yield  a  solid  angle  of  1.4963  mst.  at  the 
target. 

A. 4. 2  The  Wire  Chamber  Efficiency 

The  wire  chamber  efficiency  remained  relatively  constant  at 
around  71%  throughout  the  experiment.  If  CMWOK  stands  for  the  wire 
chamber  counts  excluding  the  misses  and  multiples  and  if  CPROT 
stands  for  the  total  counts,  the  wire  chamber  efficiency  is  given 
by: 


CMWOK 

CPROT 


A.  3 


A. 4. 3  The  System  Live-Time 

L  represents  system  live-time  and  was  calculated  using: 

Busy  Gates  -  Fast  clears  to  CAMAC  , 

L  MRS  AND  GATE  OUTPUT 

For  runs  1-8,  the  focal  plane  scintillator  was  not  in  trigger. 
From  run  9  onwards,  it  was  put  in  trigger  in  order  to  limit  the 
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value  of  This  essentially  simulated  the  right  hand  side 

electronics  and  so,  we  have  another  expression  for  system  live¬ 
time. 

t  _  Busy  Gates  -  Fast  clears  to  CAMAC  . 

RHS  AND  MRS  GATE  OUTPUT  * 

The  live-time  varies  from  14.62%  for  run  8  to  94.11%  for  run 

18. 


A. 4. 4  Beam  Current 

The  total  number  of  incident  protons  was  determined  by  the 
polarimeter  situated  upstream  from  the  T1  target.  Refer  to  Chapter 
3  for  the  beam  current  evaluation. 

The  target  thickness,  T  is  given  by: 


T 


net 


P  •  ^ 


d  a 
dfl 


A. 6 


where  N  is  the  number  of  events  due  to  tritium  alone,  P  the 
net 

number  of  incident  protons,  the  differential  cross  section  for  PT 
elastic  scattering  taken  from  Ref.  Bli-82,  S  the  solid  angle,  L  the 
system  live-time,  and  M  the  wire  chamber  efficiency. 


net 


=  N 


TiT 


A. 7 
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where : 


N  *  P 
Ti  TIT 

L  *  M  *  P 
Ti  Ti  Ti 


A. 8 


where  is  the  number  of  counts  corresponding  to  the  tritium 

region  in  the  background  run,  P  .j,  the  number  of  incident  protons 
for  the  tritium  run,  LTi  the  system  live-time  for  the  background 
run,  the  wire  chamber  efficiency  for  the  background  run  and  PTi 

the  number  of  incident  protons  for  the  background  run. 

The  results  are  compiled  in  table  A. 2.  The  average  thickness 
was  evaluated  with  an  uncertainty  of  7%. 

Tritium  content  in  the  target  =  1.220  ±  0.084  mg/cm2. 

A. 5  CORRECTION  FOR  THE  RADIOACTIVE  DECAY 

The  experiment  P  +  T  -*  a  +  y  was  done  one  year  previous  to  the 
time  when  target  thickness  evaluation  was  done.  This  calls  for  a 
correction  in  the  thickness  to  account  for  the  radioactive  decay  of 
tritium.  Tritium  undergoes  8  decay  with  a  half-life  of  12.9  years. 
If  Nq  stands  for  the  tritium  content  at  the  time  the  experiment  was 
done ,  we  have : 


N  0  =  (1.220  ±  0.084)  e  °*  693/12-  9 


0  =  1.287  ±  0.088  mg/cm2. 


i.e.  N 


TABLE  A. 2 


THE  TRITIUM  CONTENT  IN  THE  TiT  TARGET 


Scattering  Angle 
;  (°) 

Tritium  Thickness 
(mg/cm2) 

Error  in  Thickness 
(mg /era2) 

18.04° 

1.314 

0.119 

I  22.01° 

1.114 

0.124 

25.99° 

1.260 

0.382 

The  average  thickness  =  1.22  ±  0.084  mg/cm2. 


APPENDIX  B 


6-RAYS 


The  head-on  collision  between  a  proton  and  an  electron  causes 

the  electron  to  be  ejected  out  of  the  atom,  whereby  the  electron 
gets  into  the  front  detectors,  simulating  good  events.  These 
energetic  electrons  are  known  as  6-rays. 

B.l  ENERGY  OF  6-RAYS 

A  close  collision  between  a  proton  and  an  atomic  electron  may 
be  approximated  as  the  collision  between  the  proton  and  free 
electron.  Fig.  B.l  gives  the  lab  frame  situation  of  this 
collision. 

Let  m  and  m^  be  the  masses  of  the  proton  and  electron 
respectively.  P  and  P"  represent  the  momentum  of  the  proton  before 
and  after  the  collision.  P*  and  E'  are  the  momentum  and  kinetic 
energy  of  the  electron  after  collision.  P  and  P*  make  an  angle  0. 
Conservation  of  momentum  and  total  energy  give  the  following 
equation  for  the  electron  kinetic  energy: 


E' 


2m  c  2 
e 


_ _ P2C  2  cos  20 _ 

[nigC  2  +  (P  2c  2  +  m^  4)  1/2]  -  p  ^  2  cos  20 


(B. 1.1) 


A  head-on  collision  means  0  =  0  in  which  case  the  electron 

energy  is  maximum,  E*. 

m 
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E ' 
m 


2m  c  2 
e 


[n^c  2  + 


P^2 _ 

(P2c2  +  m^14) 


p2c2 


(B. 1.2) 


The  front  detectors  in  the  experiment  were  kept  at  an  angle  of 
~15°  for  300  MeV  incident  protons,  ~16°  for  227  MeV  protons  and  ~18° 
for  373  MeV  protons.  This  means  that  the  6-rays  that  had  to  be 
tolerated  in  the  absence  of  the  magnet  were  of  energies  0.7  MeV,  0.5 
MeV  and  0.8  MeV  respectively. 

Since  m  »  ,  eq.  (B.1.1)  reduced  to  the  following  form: 

2m  c  2 

E’  =  - - ~ -  (B.1.3) 

1 

P 

where  8  =  t; -  ,  E  being  the  total  energy  of  the  proton.  We  also 

L  tOt 

tot 

have : 


2m  c  282 

E. - ~  —  (B.1.4) 

m  1  -  32 

B. 2  6-PRODUCT ION  PROBABILITY 

The  generalization  of  the  Rutherford  formula  to  account  for  the 
fact  that  the  electrons  produced  are  relativistic,  yields  the  cross 
section  for  6-production.  The  Rutherford  formula  takes  the  force  of 
interaction  from  the  impact  parameter,  finds  the  momentum 
transferred  from  the  force  and  the  collision  time,  and  then  gets  the 
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kinetic  energy  transferred  to  an  electron  from  the  momentum, 
assuming  the  non-relat ivistic  case.  Then  it  makes  use  of  the 
probability  of  collision  with  the  impact  parameter  being  between  b 
and  b  +  db  in  the  thickness  dx  gm  cm'2  as  follows: 


F(b)db  dx  =  2irb  db  N  ^  dx 


(B.2.1) 


where  N  is  the  Avogrado  number,  Z  the  atomic  number  and  A  the  atomic 
mass  number. 

Since  the  energy  transferred,  E',  is  already  obtained  in  terms 
of  impact  parameter,  probability  can  be  expressed  in  terms  of  E'. 


2B  m  c2 


p  HE' 

$  ( E  *  )  dE  *  dx  = - - - - -  dx 

32  (E')2 


(B.2.2) 


This  is  the  probability  for  the  6-ray  of  energy  E'  to  be 
produced  while  the  proton  travels  in  the  material  a  distance  dx 
gm  cm”2.  The  factor  B  represents  the  total  area  covered  by  the 
electrons  contained  in  one  gm,  each  electron  considered  a  sphere  of 
radius  r  . 


e 


Z  o 
B  =  ttN  — r —  r  2 
A  e 


(B.2.3) 


where  r  = 
e 


,  the  classical  electron  radius,  which  equals 


2.8  fm. 
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For  the  relativistic  case,  which  is  the  situation  in  our  case, 
eq .  (B.2.2)  generalizes  to: 


$  (E,E ’ )  dE» 


2Bm  c 
e 


dE 


(E') 


[i 


E+mc 


;)  ]  (B.2.4) 


where  E'  and  E^  are  given  by  eqs.  B.1.3  and  B.1.4.  Here  B  signifies 
the  collision  probability  dependence  on  the  electron  density  and 
gives  the  dependence  of  the  energy  transfer  on  the  collision  time. 

Eq.  B.2.4  may  be  rewritten  in  terms  of  the  scattering  angle  0. 
We  get: 


$  (E ,  0)  d  0  = 


2B  sec  20  tan  0 
34  (sec  20-  32) 


-  (sec 20  -  2 32  +  34)  d0 


(B.2.5) 


Integrating  this  over  the  angles  spanned  by  the  wire  plane,  it 
is  found  that  the  expected  6-rays  flux  for  300  MeV  protons  incident 
on  the  TiT  target  would  be  5.6  x  105/s. 


APPENDIX  C 


CHERENKOV  COUNTER  CALIBRATION  AND  EFFICIENCY  DETERMINATION 


The  purpose  of  this  experiment  was  to  determine  the  efficiency 

of  the  Cherenkov  counter  for  detecting  positrons.  The  response  of 
the  counter  for  positrons  of  various  energies  was  determined.  These 
data  on  'e'  detection  efficiency  were  then  used  in  calculating  the 
efficiency  for  detecting  y  rays,  knowing  the  correspondence  between 
the  y  energy  and  the  e+  or  e  energy  from  pair  production  and 
Compton  scattering.  This  calculation  of  y  detection  efficiency 
assumed  that  the  'ef  detection  efficiency  is  independent  of  its 
position  in  the  detector.  However,  further  calculations  were  then 
made  to  take  into  account  the  shower  loss,  e  loss,  or  both  from  the 
knowledge  of  the  detector  geometry  as  well  as  the  knowledge  of  the 
detector  positioning  with  respect  to  the  target.  The  distribution 
of  the  energies  of  e+  and  e  for  a  particular  y  was  determined  using 
a  Monte-Carlo  simulation  program. 

C.l  M20  STOPPED  tt/  p  CHANNEL 

The  positron  beam  available  at  the  M20  stopped  tt/  jj  channel  of 
the  TRIUMF  Cyclotron  Facility  was  used  for  the  experiment.  M20  is  a 
9  m  long  secondary  channel  available  at  the  target  station  1AT2. 

The  TRIUMF  Facility  has  two  pion  production  target  stations,  1AT1 
and  1AT2,  which  are  situated  in  beam  line  1A.  The  target  at  1AT2 
was  a  5  cm  Be  when  the  experiment  was  done.  The  tt+,  generated  at 
the  target,  decays  at  rest  in  the  surface  layer  of  the  target,  in 
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the  mode  tt+  p+u.  The  produced  p+  underwent  decay  in  the  mode  p+  + 
e+uu.  Typical  values  of  the  fluxes  of  ir+,  y+  and  e+  for  a  1  pA 
proton  beam  of  500  MeV  at  1AT2  are  106/s,  4  x  10  Vs  and  5  x  105/s 
respectively.  A  set  of  dipoles  and  quadupoles  provided  momentum 
selection  up  to  17.5  MeV/c  with  a  resolution  of  5%.  The  Cherenkov 
counter  detected  tt+,  p+  and  e+,  but  e+  events  were  separated  from 
the  rest  of  the  events  by  the  TOF  measured  between  the  RF  of  the 
cyclotron  and  a  counter  in  the  beam.  RF  provided  the  timing 
required. 

C.2  THE  EXPERIMENTAL  DETAILS 

The  purpose  of  the  experiment  was  to  find  the  efficiency  of  the 
Cherenkov  counter  for  detecting  e+  of  momenta  varying  from  12  MeV/c 
to  175  MeV/c.  The  response  of  the  counter  for  e  beam  incident  at 
various  positions  across  the  face  of  the  detector  was  also  measured. 
The  latter  measurement  was  necessary  to  get  an  idea  of  the  shower 
loss  at  the  edges  of  the  counter. 

Fig.  C.l  gives  the  electronics  used  in  the  experiment.  C 
represents  the  Cherenkov  counter,  5"  diameter  x  7"  long  (corresponds 
to  7.5  radiadion  lengths  for  lead  glass),  P  is  the  plastic 
scintillator,  6"  x  6",  T  is  the  tagging  plastic  scintillator,  1"  x 
1" .  Looking  along  the  direction  of  the  beam,  C  was  preceded  by  P 
which  was  preceded  by  T.  The  distance  between  C  and  P  was  1  1/2", 
that  between  P  and  T  was  3”.  The  level  of  the  discriminator 
attached  to  C  was  adjusted  using  an  attenuator  in  the  anode  circuit 
such  that  at  17  db  attenuation,  the  count  rate  from  an  Americium 
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source  noted  by  a  rate  meter  (ORTEC  441)  was  half  of  the  count  rate 
at  0  db.  This  discriminator  level  was  the  one  present  during  the 
PT  ay  experiment.  The  ratio  of  SCALER  3  to  SCALER  2  represented 
the  beam  focussing  quality.  The  ratio  of  SCALER  1  and  SCALER  4 
represented  the  Cherenkov  counter  efficiency. 

The  experiment  was  performed  at  electron  momenta  12,  15,  25, 

37,  50,  100  and  175  MeV/c.  The  measurements  were  made  with  the  beam 
incident  at  distances  0,  2,  4,  6  and  8  cm  from  the  center  of  the 
Cherenkov  counter  face. 

C.3  THE  RESULTS 

The  central  efficiency  was  13%  at  12  MeV/c  (e+  momentum)  which 
increased  rapidly  to  56%  at  15  MeV/c,  95%  at  25  MeV/c  and  remained 
at  100%  at  and  beyond  50  MeV/c.  Positional  response  was  such  that 
for  a  100  MeV/c  e  beam  incident  at  8  cm  from  the  center,  the 
efficiency  was  45%;  for  6  cm,  75%;  for  4  cm,  91%;  and  for  1  cm, 

100%.  Fall-off  of  the  efficiency  with  respect  to  distance  from  the 
center  was  more  rapid  for  decreasing  energy  (Hut-80).  An  empirical 
expression  for  the  Cherenkov  counter  efficiency  as  a  function  of  e 
momentum  was  found  to  be  (Ste-80) 

E  =  1  -  e"™  (C.3.1) 

where  m  =  0.591  . 

E  is  the  efficiency;  p  the  e  momentum  in  MeV/c. 


' 
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C.4  MONTE-CARLO  CALCULATION  FOR  PHOTON  DETECTION  EFFICIENCY 

The  photon  entering  the  lead  glass  triggers  a  long  series  of 
interactions  resulting  in  a  flux  of  electrons,  positrons  and 
photons,  known  as  an  electromagnetic  cascade  shower.  It  is  easier 
to  explain  this  phenomenon  through  a  simulation  program  than 
mathematically,  hence  the  need  for  a  Monte-Carlo  calculation. 

Photon  interactions  with  matter  may  be  classified  into  three 
categories.  The  interaction  of  a  photon  with  an  atom  as  a  whole 
leads  to  a  photoelectric  effect  where  the  electron  is  knocked  out  of 
the  atom.  This  means  that  the  photon  energy  has  to  be  at  least  the 
same  as  the  electron  binding  energy.  The  interaction  of  the  photon 
with  a  free  electron  leads  to  the  Compton  effect  where  most  of  the 
photon  energy  is  handed  over  to  the  electron.  The  photon  energy  is 
reduced  by  the  kinetic  energy  of  the  recoil  electron.  The 
interaction  of  the  photon  with  the  coulomb  field  of  the  nucleus  is 
known  as  pair  production  where  the  photon  disappears  giving  rise  to 
an  electron  and  positron.  Obviously,  the  threshold  condition  is 
that  the  photon  energy  should  exceed  the  sum  of  rest  energies  of  the 
electron  and  positron.  The  excess  energy  appears  as  kinetic  energy 
of  the  electron  and  positron.  The  recoil  nucleus  accounts  for  the 
momentum  balance.  At  y  energies  of  interest  to  us  (100’s  of  MeV), 
only  the  Compton  effect  and  pair  production  are  important.  Again, 
as  the  energy  increases,  the  Compton  effect  becomes  less  important 
and  pair  production  becomes  more  dominant. 

The  charged  particles,  ie.  e+  and  e“  produced  in  the  lead  glass 
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by  pair  production  as  well  as  Compton  scattered  are 
ultrarelativistic.  That  is,  they  travel  in  the  medium  with  a 
velocity  greater  than  the  velocity  of  light  in  the  medium.  This 
generates  radiation  along  a  conical  wavefront  known  as  Cherenkov 
radiation.  The  angle  of  emission  of  the  radiation,  0,  is  given  in 
terms  of  the  refractive  index  of  the  medium  'n'  and  *  3*  of  the 
particle  in  the  following  way: 


cos 


(C.4.1) 


For  the  lead  glass  which  constituted  the  material  of  the 
Cherenkov  counters,  n  =  1.67.  The  value  of  M  1  which  gives 
0  »  53°. 

The  number  of  Cherenkov  photons,  N,  emitted  per  unit  length  of 

dN 

the  radiator  is  denoted  by  -p-  . 

y  dx 


dN  =  2  nog  2  (i _ 1 _ 

dx  c  n232 


(C.4.2) 


where  a  is  the  fine  structure  constant,  1/137,  Z  is  the  atomic 
number  of  the  material  through  which  the  charged  particle 
traverses. 

The  charged  particle  generated  in  the  pair  production  process 
gets  deflected  in  the  coulomb  field  of  the  nucleus  and  hence 
experiences  acceleration  77,  where  Z  is  the  atomic  number  of  the 
nucleus  and  M  the  electron  mass.  According  to  classical 
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electrodynamics,  the  resulting  radiation,  known  as  brerasstrahlung 

Z  2 

radiation,  has  an  intensity  proportional  to  —  .  An  electron  of 

M2 

kinetic  energy  T  generates  a  photon  of  any  energy  up  to  T. 

In  general,  we  have  the  following  situation  triggered  by  an 
incoming  energetic  photon  in  lead  glass.  It  may  either  pair  produce 
or  Compton  scatter.  If  pair  production  takes  place,  the  generated 
e+  and  e”  produce  bremsstrahlung  photons  which  can  start  a  new  chain 
of  interaction.  If  Compton  scattering  takes  place,  (i)  the  recoil 
electron  will  produce  a  photon  by  bremsstrahlung  and  (ii)  the 
Compton  scattered  photon,  if  energetic  enough,  will  start  a  new 
chain  of  interaction.  The  photon  generated  in  process  (i)  can  also 
trigger  yet  another  chain  of  interactions.  This  results  in  an 
electromagnetic  cascade  shower,  the  flux  depending  upon  the 
radiation  length  traversed.  The  Cherenkov  light  produced  is 
proportional  to  the  total  path  traversed  by  both  e+  and  e“,  known  as 
track  length,  which  in  turn  is  proportional  to  the  energy  of  the 
primary  photon.  This  means  that  the  total  Cherenkov  light  collected 
is  proportional  to  the  energy  of  the  primary  photon,  provided,  of 
course,  that  the  entire  shower  is  contained  in  the  counter 
material. 

A  typical  cascade  simulation  program  will  concern  itself  with 
things  like  the  distribution  of  energies  in  the  pair  production,  the 
distribution  of  energies  of  photon  and  electrons  in  the  Compton 
scattering  and  brerasstrahlung,  and  assign  these  values  arbitrarily 
from  the  interaction  probabilities  known  from  the  energies  of  the 
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particles  involved  and  the  mass  number  and  atomic  number  of  the 
radiator  material  elements  (Lon-75).  Moreover,  it  will  take  into 
account  the  shower  loss  which  is  a  function  of  the  geometry,  the 
term  geometry  standing  for  the  shape  and  size  of  the  radiator 
material  as  well  as  the  positioning  of  the  counter  with  respect  to 
the  target.  The  shower  fluctuations  from  track  length  fluctuations 
generated  by  the  geometry  could  cause  poor  resolution  in  the  photon 
energy  detected.  Since  at  large  shower  depths,  there  are  relatively 
large  numbers  of  low  energy  electron,  collision  loss  is  more 
dominant  than  radiation  loss.  Multiple  coulomb  scattering  of  the 
electrons  by  atoms  as  well  as  the  Compton  process  gives  rise  to 
lateral  spread  of  the  shower.  The  simulation  program  could  also  try 
to  give  an  estimate  of  the  lateral  spreading. 

The  Monte  Carlo  simulation  program,  written  to  calculate  the 
photon  efficiency  of  the  counters  used  in  the  experiment  PT  -*■  ay,  is 
described  in  (Ste-80).  This  program  is  in  the  spirit  of  the 
calculations  done  by  D.I.  Sober  (Sob-73).  A  brief  explanation  of 
(Sob-73)  is  warranted. 

Sober  takes  pair  production  and  Compton  scattering  as  the 
interactions  which  the  photon  undergoes.  The  probabilities  are 
summed  over  the  lead,  silicon  and  oxygen  atoms  in  the  lead  glass. 

The  detection  probability  of  a  photon  of  energy  k  having  traversed  a 
distance  t  is  expressed  as  the  product  of  probability  of  interaction 
and  probability  of  detecting  one  or  more  'ef  produced  in  the 


interaction: 
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ey  (k,t)  -  (1  -  e  u(k)t)  E  (k,t)  (C.4.3) 

where  the  first  factor  is  the  probability  that  the  photon  has 
interacted  in  the  counter  before  reaching  * t  * .  The  second  factor  is 
the  average  efficiency  for  detecting  one  or  more  electrons.  This 
factor  was  obtained  from  the  Monte-Carlo  calculation  with  the 
following  assumptions:  (i)  the  electron  detection  efficiency  is 
independent  of  its  position  within  the  counter,  (ii)  the  position  of 
interaction  of  the  photon  and  distribution  of  final-state  energies 
are  assigned  by  the  Monte-Carlo  technique,  (iii)  the  Compton 
scattered  photon  is  allowed  to  propagate  further  provided  the  energy 
is  greater  than  2  MeV,  and  (iv)  after  a  pair  event,  the  shower  is 
not  propagated  further.  The  efficiency  of  detecting  the  electron 
was  taken  from  the  electron  calibration  data. 

The  calculations  described  in  (Ste-80)  took  the  empirical 
expression  C.3.1  for  the  electron  detection  efficiency.  Again  only 
Compton  scattering  and  pair  production  were  considered.  An  event 
was  terminated  at  pair  production.  Compton  scattering  was  followed 
up.  The  event  was  terminated  if  the  secondary  photon  pair  produced 
left  the  detector  or  attained  the  subthreshold  energies. 

As  may  be  seen,  the  loss  of  the  electrons  and  positrons  were 
still  to  be  accounted  for.  This  was  done  in  (Hu-80)  where  the  track 
length  traversed  was  calculated  as  a  function  of  the  solid  angle 
subtended  at  the  target  by  various  annular  regions  of  the  Cherenkov 
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counter  face.  The  radiation  length  to  be  traversed  by  the  electron 
in  order  to  produce  enough  light  to  be  seen  by  the  photomultiplier 
was  calculated,  taking  into  account  the  discriminator  level,  and  the 
loss  of  electrons  at  the  edges  of  the  detector  was  estimated. 

The  results  of  the  calculations  after  taking  into  account  both 
the  photon  losses  and  electron  losses  are  tabulated  in  Table  C.l. 
Values  of  the  efficiency  for  different  y  energies  and  different 
detector  distances  are  plotted  in  fig.  C.2. 
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TABLE  C.l 

PHOTON  DETECTION  EFFICIENCY  OF  LEAD-GLASS  CHERENKOV  COUNTER 


E 

Y 

EFFICIENCY  (%)  AT  THE  DISTANCES  (cm) 

(MeV) 

30 

75 

100 

100 

81 

85 

89 

200 

83 

87 

91 

300 

83 

89 

93 

400 

86 

90 

94 

204 


0.95 


205 


O 


APPENDIX  D 


STUDY  OF  THE  REACTION  p  +  d  ♦  %e  +  y 

The  motivation  for  looking  at  this  reaction  is  described  in 
section  2.2.  The  reaction  had  been  studied  thoroughly  earlier 
(Cam-82,  Cam-83)  and  therefore  provided  a  means  to  check  the 
functioning  of  the  experimental  set  up  by  measuring  the  cross 
section  and  comparing  the  result  against  the  previously  measured 
value.  Furthermore,  as  mentioned  in  section  3.4.3,  it  provided  a 
way  to  calculate  the  Cherenkov  counter  efficiency  as  a  function  of 
the  pulse  height  cut,  by  providing  a  well-defined  Cherenkov  pulse 
height  spectrum  with  good  statistics. 

D. 1  THE  EXPERIMENTAL  SET-UP 

The  same  experimental  set-up  was  used  as  for  P  +  T  a  +  y 

(Chapter  2).  A  deuterated  polyethylene  (CD^)  target  of  73.3  mg/cm2 

thickness  was  used.  The  beam  energy  was  maintained  at  3U0  MeV.  The 

recoil  particle  (3He)  detector  assembly  was  situated  at  14°  on 

either  side  of  the  beamline.  The  distance  was  43"  from  the  target. 

This  detector  assembly  identified  3He's  of  energies  75  MeV,  100  MeV 

and  125  MeV.  The  y-detectors  (veto-Cherenkov  counter  assembly)  were 

kept  at  angles  49°,  75°  and  107°  at  distances  of  24",  20"  and  20" 

respectively  from  the  target  on  either  side  of  the  beam  line.  These 

viewed  y’s  of  energies  224  MeV,  199  MeV  and  173  MeV.  The  above 

CM 

angles  corresponded  to  0^  =  62°,  90°  and  120°.  Distances  were 

adjusted  such  that  the  Cherenkov  counter  image  fell  totally  within 
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the  wire  planes,  which  meant  that  the  solid  angle  was  determined  by 
the  Cherenkov  counter.  The  total  collected  charge  was  36  jjc.  The 
differential  cross  section  and  the  analyzing  power  were  measured  at 
the  above  three  points. 

D. 2  THE  DATA  ANALYSIS 

Data  for  this  run  did  not  need  to  be  rewritten  because  there 
was  no  problem  of  3He  stopping  in  the  front  plastic  and  so  LEFT  and 
RIGHT  OCR  tests  were  unambiguous.  The  analysis  procedure  was  the 
same  as  described  in  Chapter  3.  The  Cherenkov  efficiencies  were 
about  88%,  the  wire  chamber  efficiency  72%,  and  the  system  live-time 
76%.  Fig.  D.l  gives  the  plot  of  recoil  particle  energy  versus 
recoil  angle.  The  distinct  peak  due  to  (p,y)  may  be  seen.  The 
carbon  subtraction  was  done  by  describing  a  box  (of  the  same  size  as 
for  (p,y)  events)  in  the  general  background  region  of  the  spectrum 
and  subtracting  those  counts  from  (p,y)  counts.  A  typical  Cherenkov 
spectrum  is  shown  in  fig.  3.8. 

D.3  THE  RESULTS 

The  cross  section  and  the  analyzing  power  were  determined  for 

r»  u 

£r  =  62°  and  90°.  Part  of  the  Cherenkov  image  for  the  angle 
CM 

0^  =  120°  fell  outside  the  wire  chamber  frame  and  that  point  was 

omitted.  However,  analyzing  powers  were  calculated  for  0^  =  120°. 
Table  D.l  gives  the  run  summary.  The  cross  section  and  the  standard 
deviation  were  calculated  using  equations  3. 3. 9. 3  and  3.6.1,  where 
the  variables  were  replaced  by  those  for  the  CD2  run.  The  analyzing 
power  and  the  error  were  evaluated  using  equations  3.7.1  and  3.7.3. 


208 


(p,Y)  events 


209 


—  —  —  cm  m 


co 

to  cj 

i  —  N 

'  cm 


CM  <*i 
—  fs!  — •  — 


(SI 


CO  —  —  CO  —  CM  CO  — 

— >  —«  <H  -s  —  —  — 

—  —  —  cJ  H  C*l  -  -  CM  — 

—  ~  rV  —  m  —  — ■  cm  — 

NS  St  CM  —  —  —  —  — 

cm  uo  co  c»  ■*.  fv'  —  —  —<  cm  —  — i 

r*  CM  —  CM  nNvC'lUMO-iN-"H  — < 

—  —  —  CMCO'<r<;cocr.  rncMLn  -  —  —  cm  —  —  —  — < 

—  —  co  — >  in  NTH  CJ  — ■  O  <  — <  CM  —  —  CM 

— >  CO  N'  L-J  li!  _J  +  CL  I  m  r  st  »-  CO  — I  CM  CO  — 

cm—*  cm  n  xi  u.  c.  I- W  _i  X  in  o  «  co  —  _  —  —  cm  cm 

—  CM  —  w-L'JlLO.Oh  +  -IU)sf  -i  CM  CON--  CM 

—•  -CMCOOO  +  +  QLJlO-C\l  CM  —  —i  CM  — 

— i  «-<  CO  NCOC3UO  +  +  QClJJCJ'C-'  —  — 

—  —  CM  CM  —  —  CO  +  +  +  +  r  —  — 

_  —  —i  ^i^fL"os++mr^  co  —  —  —  —  — >  — 

'  CM  — <  C0C0C0'OO=>r0>mCM  CsJ 

I  — ■  CM  "C  — i  Csl  —  i/»CJLi_1jJLOCQCsICM  —  —  — < 

—  — —  -".<iiuD<n--  — 

—  —  CO  CM  —  —  COCOG~i<^SDLfi  —  cm  —  —  — 

— «  —  —  CMCMLniosomrsin  —  — 

1  —  —  '<t  —  — •  in  co  in  ■m-  co  —  —  — «  — 

—  CM  —  —  —  CM  —  CMUO'#''*  —  CM  —  —  —  —  — 

—  — >  CM— <  CM  —  —  CM  —  Csl  CM  —  —  — <  —I 

— <  — •  —  Cv!  —  ——————  t\|  —  CM 

— <  — —  — <  — <  CO  CM  — <  0">  CM  —  CO  — <  — < 

CM  —  CM  —  — >  "<  —  Csl  —  —  — > 

CM  —  —  — «  CO  —  CM  —  —  — ■  —  —  — <  CM  — • 

— <  —  —  CM  CO  —  CM  —  —  CM  — ■  CO  Csl  V  CO  CM  — t  CM  CM  — 

\  —  —  «<•  CO  — •  — <  -•  CM  CO  —  —  —  — <  — I 

—  —  — <  —  —  —  CM  CM  CM  — <  CM  CM 

—  — <  —  —  —  CM  —  —  —  —  —  —  CM 

— t  CM  CO  CO  CM  —  —  —  CM  —  CM  CM 

—  —  CM  — •  — •  CM  —  C>  —  —CM  — *  —  — 

—  — t  CM  — .  — •  — i  CO  —  — «  CO  CM  •«(  CM  — •  — < 

CM  — i  CM  CM  —  —  —  CM  CO  — ,  —  — ,  —  — 

<*>—„._  —  CM  —  — <  CO  — «  — > 

CM  CM  —  —  —  —  -CO  —  — <  CO  CM  — >  CM  —  CM 

— <  CM  — i  — •  — <  —  — <  —  —  — •  —  —  — 

—  CM  —  CM  —  CM  —  —  —  CM  CO 


CM  CM  —  —  CO  —  — 
Cs!  —  —  CM  Csl  —  —  — 
CO  CM  CM  CO  — 

—  ri  C)  SO  -C  CM  ■<  CO  — 
LH  CO  CM  —  CO  —  CM 


—  —  CM 


CO 


A0tf3N3  aiDITNVd  TE0D3N 


bj 

b< 

♦ 

w 

•  ••  •  . . fc<  tsi  bj  ti  t?  tx.  t-g  to  Isa  tv 

OOb  OD^OOOOOCHsOtsitsitw&tsttsCjCjisiislsiisilsOfsits . tsj 

Ov%ftN  o  ic  -m  r*»oi»*cj<T\corsi£>uv^c*>  cm  —  t?cncorMvmn'4  co  cm  —  &jfe>lsC5t^t£)<abjCsiCS)Csj 

roror0Cfi»*0COrOfnr-»co(MCMrj(MrMCMiCMCMCsiftl  —  —  —  —  —  —  —  —  —  —  cr.cor^-uOLn>*ccofM—  • 

b5 


RECOIL  PARTICLE  ANGLE 


210 


TABLE  D.l 

P  +  d  +  %e  +  y  DATA  SUMMARY,  (P,y)  EVENTS 


LEFT 


Ll  (e™=  62°) 

L2  (e^M  =  90°) 

SPIN  UP 

SPIN  DOWN 

SPIN  OFF 

SPIN  UP 

SPIN  DOWN 

SPIN  OFF 

697  ±  26 

792  ±  28 

312  ±  18 

634  ±  25 

671  ±  26 

263  ±  16 

RIGHT 


R1  (0™=  62°) 

R2  (9CM  =  90°) 

y 

SPIN  UP 

SPIN  DOWN 

SPIN  OFF 

SPIN  UP 

SPIN  DOWN 

SPIN  OFF 

779  ±  28 

689  ±  26 

270  ±  16 

576  ±  24 

581  ±  24 

213  ±  15 

Table  D.2  summarizes  results  obtained  after  the  accidentals 
subtraction.  Fig.  D.2  plots  the  cross  section  values  along  with  the 
existing  values  (Cam-83).  Fig.  D.3  gives  the  plot  of  the  analyzing 
power.  The  results  (fig.  D.2)  may  be  seen  to  be  consistent  with  the 
previously  existing  data. 
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TABLE  D. 2 

CROSS  SECTION  AND  ANALYZING  POWER  FOR  THE  REACTION  P  +  d  -►  3He  +  y 

AT  T  =  300  MeV. 

E 


5  >- 

CD 

Cross  Section 
(nb/sr ) 

Analyzing  Power 

62° 

61.2  ±  1.3 

0.08  ±  0.07 

90° 

48.0  ±  1.0 

0.01  ±  0.10 

120° 

— 

-0.03  ±  0.16 

FROM  To SE. 


